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ABSTRACT 


From  a  customers  point  of  view,  warranties  seem  desirable  as  a 
means  of  improving  reliability  (availability) and  reducing  overall  system  cost. 
However  contractors  may  be  reluctant  to  provide  warranties,  particularly 
when  an  assumption  of  excessive  risk  is  perceived  to  exist.  In  such 
situations  a  mutually  agreeable  arrangement  may  be  achieved  by  means  of 
an  incentive  fee  type  of  warranty.  This  point  of  view  is  similar  to  the 
various  incentive  contract  forms  introduced  some  twenty  years  ago  to  divide 
risk  assumption  between  the  customer  and  the  contractor.  In  incentive 
contracts  it  is  usual  to  have  some  linear  decrease  in  fee  as  cost  increases. 

To  carry  the  incentive  concept  into  the  structuring  of  warranties  a  different 
point  of  view  is  required,  based  on  the  same  general  concepts  of  dividing 
risk  assumption.  Incentivising  reliability  warranties  is  achieved  by  pro¬ 
viding  a  warranty  fee  pool  as  a  function  of  time.  Use  of  an  increasing  fee 
pool  provides  the  incentive  for  a  contractor  to  assume  the  risk  involved  in  a 
warranty  situation  by  sharing  risk  with  the  customer.  It  is  also  possible  for 
the  contractor  to  achieve  additional  profit  under  properly  structured 
incentivised  warranties. 

^'^he  research  reported  here  defies  the  concepts  of  incentivised 
reliability  (availability)  warranties,  develops  mathematical  models  of  the 
structure  of  such  warranties,  and  analyzes  selected  examples.  The  goal  of 
the  analyses  is  to  identify  and  study  the  role  of  contract  variables  as  they 
relate  to  the  division  of  risk  assumption,  the  improvement  of  reliability 
(availability)  and  the  control  of  cost  in  a  warranty  situation. 

Since  the  research  deals  with  reliability  considerations  it  is  carried 
out  in  a  stochastic  framework,  where  the  mathematical  models  provide 
expected  values  and  standard  deviations  of  costs  for  repair/replacement  as 
functions  of  time.  Warranty  fee  pools  are  represented  by  various  functional 
forms  and  different  probability  laws  are  employed  to  govern  the  failure  events 
and  times  to  failure.  / 
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From  the  point  of  view  of  the  user,  delivery  of  a  piece  of  equipment 
is  only  the  beginning  oi  a  period  over  which  the  equipment  is  required  to 
function,  its  lifetime.  On  the  other  hand  a  contractor  or  supplier  may  view 
the  delivery  as  the  end  of  its  period  of  concern  with  the  equipment.  These 
two  viewpoints  do  not  represent  desirable  approaches  to  equipment  acquisi¬ 
tion  and  in  actual  practice  the  user  and  supplier  broaden  their  periods  of 
interest  in  the  equipment,  both  being  involved  throughout  at  least  part  of  the 
development  and  the  use  periods.  User  involvement  in  development  occurs 
through  design  specifications,  design  reviews,  and  various  acceptance  test 
procedures.  Supplier  involvement  in  equipment  operation  occurs  in  some 
form  of  maintenance  agreement  to  be  found  in  most  contracts  for  equipment 
acqusition.  One  way  to  achieve  satisfactory  use  is  by  improving  the  reli¬ 
ability  of  equipment.  However  reliability  is  a  difficult  quantity  to  design, 
being  composed  of  various  random  features.  A  supplier  often  agrees  to 
some  level  of  reliability  as  part  of  a  contract  but  in  fact  the  actual  level  of 
reliability  may  be  found  to  differ  from  specified  values  once  the  equipment 
is  operational.  In  addition,  being  random,  the  times  to  failure  of  equipment 
may  prove  to  be  considerably  short  while  still  falling  technically  within 
design  specifications. 

In  recent  years  there  has  been  an  attempt  to  improve  reliability  by 
having  the  contractor  assume  the  risk  of  an  equipment  failure  in  terms  of  a 
warranty  on  the  equipment.  The  Reliability  Improvement  Warranty  (RIW) 
concept  has  received  considerable  attention  and  analysis,  as  typified  by  the 
studies  reported  in  (1).  One  of  the  main  features  of  an  RIW  is  that  it  forces 
a  contractor  to  take  a  serious  view  regarding  reliability.  This  is  not  done 
in  many  cases  because  of  the  difficulty  associated  with  quantifying  reliability 


in  a  significant  way.  Data  on  reliability  are  apt  to  be  scarce  or  inaccurate. 
The  concepts,  based  on  probability  theory,  seem  strange  to  most  manage¬ 
ment  and  engineering  personnel,  and  therefore  only  cursory  attention  may 
be  paid  to  contractual  specifications  in  the  reliability  area.  The  RIW  makes 
reliability  considerations  important  from  a  monetary  point  of  view,  acting 
as  a  form  of  penalty  against  profits  should  the  equipment  need  to  be  re¬ 
paired  or  replaced  during  the  warranty  period. 

From  the  users  point  of  view  the  RIW  has  two  important  features. 

It  focuses  the  contractor  on  reliability  considerations  which  are  likely  to 
improve  the  operational  status  of  equipment,  and  requires  the  contractor  to 
share  the  risk  of  failing  to  achieve  a  satisfactory  level  of  operation. 

The  RIW  concept  introduces  two  possibilities  for  generalization. 
Though  it  is  desirable  from  the  users  point  of  view,  contractors  may  be 
very  reluctant  to  assume  the  full  level  of  risk  implied  in  a  normal  warranty 
contract.  This  leads  to  the  consideration  of  various  incentive  forms  of 
warranties  in  which  the  risk  can  be  shared  between  the  contractor  and  the 
user.  Furthermore  in  such  situations  there  can  exist  the  possibility  of 
additional  profit  for  the  contractor  when  high  reliability  is  achieved.  Since 
such  achievement  is  often  the  major  goal  of  the  user  as  well,  he  has  a  strong 
motivational  instrument  in  the  properly  structured  incentive  warranty.  The 
other  direction  of  generalization  is  to  the  concept  of  availability  which,  in 
fact,  is  most  often  the  users  concern  rather  than  purely  reliability  alone  * 
Equipment  failure  can  often  be  tolerated  so  long  as  down  time  is  sufficiently 
short  and  the  time  between  failures  satisfactorily  long. 

This  report  studies  the  concept  of  incentivisation  of  warranties  for 
die  reliability  and  more  generally,  availability  situations.  The  methodology 
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is  based  on  mathematical  models  of  warranty  contract  structures  in  a 
stochastic  framework  as  described  in  this  chapter. 

Section  1.  1  Introduction 

Incentive  type  contracts  have  been  in  use  for  some  time.  They  are 
designed  to  motivate  a  contractor  to  keep  cost  low,  meet  schedules,  or 
achieve  high  levels  of  performance.  The  most  common  form  of  incentive 
contracts  provide  for  a  fee  that  decreases  with  increasing  cost.  In  some 
cases  the  decrease  may  reach  to  a  negative  fee  (penalty),  whereas  other 
forms  of  contract  limit  the  decrease.  Most  often  the  decrease  in  fee  is 
linear  or  consists  of  a  few  linear  segments  differing  in  their  rate  of  de¬ 
crease.  Incentive  contracts,  in  their  various  forms,  have  bden  studied 
extensively  and  widely  applied.  Theoretical  analysis  of  such  contract 
forms  may  be  found  in  (2)  and  (3). 

The  incentive  fee  concept  was  developed  as  an  instrument  for  causing 
the  contractor  to  assume  part  of  the  risk  in  situations  where  cost  was  dif¬ 
ficult  to  determine.  Incentivised  contracts  fall  between  fixed  fee  situations 
where  costs  should  be  deterministic  and  purely  research  type  areas  where 
cost  plus  fixed  fee  contracts  are  required  before  a  contractor  will  engage  in 
the  necessary  work.  Thus  it  is  the  relative  uncertainty  of  cost  that  makes 
an  incentive  fee  contract  desirable  to  both  contractor  and  customer.  Such 
concepts  may  be  applied  to  the  reliability /availability  situation.  If  the 
availability  profile  of  a  piece  of  equipment  was  deterministic  a  contractor 
could  figure  a  meaningful  addition  to  cost  for  including  a  warranty.  At  the 
other  extreme,  if  the  availability  behavior  was  completely  erratic  no  con¬ 
tractor  would  be  willing  to  provide  a  warranty.  Thus  the  RIW  situation  can 
only  be  agreed  to  when  a  contractor  believes  he  can  correctly  estimate  the 
availability  profile  and  include  sufficient  cost  to  compensate  his  assumption 
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of  the  warranty.  This  situation  may  not  be  expected  to  occur  too  often, 
particularly  if  the  contractor  fully  understands  the  random  character  of 
availability.  However  the  use  of  incentives  described  above  suggests  a 
method  for  greatly  enlarging  the  range  of  situations  in  which  a  contractor 
may  engage  in  a  warranty.  Rather  than  having  only  the  two  extreme  situa¬ 
tions  of:  no  warranty,  with  user  assuming  all  the  risk,  and  normal  warranty, 
with  contractor  assuming  all  the  risk,  the  incenti vised  warranty  allows  a 
distribution  of  risk  assumption  between  user  and  contractor. 

To  carry  over  the  incentive  concept  into  the  structuring  of  warranties 
a  different  point  of  view  is  required,  though  still  following  the  same  general 
principle  of  dividing  risk  assumption.  Incentivising  reliability  warranties 
is  achieved  by  providing  a  warranty  fee  pool  as  a  function  of  time.  Use  of 
an  increasing  fee  pool  provides  the  incentive  for  a  contractor  to  assume  the 
risk  involved  in  a  warranty  situation  by  sharing  risk  with  the  customer.  An 
increasing  fee  pool  motivates  the  contractor  to  increase  time  between 
failures  as  much  as  possible  which  is  one  major  feature  of  high  availability. 

It  is  also  possible  for  the  contractor  to  achieve  additional  profit  under  a 
properly  structured  incentivised  warranty.  However  such  situations  must  be 
treated  very  carefully  to  prevent  motivating  the  contractor  to  allow  failure, 
which  is  counter  to  a  major  feature  of  high  availability  i.  e.  few  failures.  A 
major  reason  for  the  detailed  mathematical  analysis  of  incentivised  avail¬ 
ability  warranties  is  to  clearly  define  the  various  features  implicit  in  such 
relatively  complex  contract  structures  and  to  quantify,  within  postulated 
stochastic  frameworks,  the  levels  and  forms  of  incenti  vis  ation  that  will 
provide  desirable  availability  profiles  at  a  cost  agreeable  to  both  the  con¬ 
tractor  and  the  customer. 
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Section  1.2  Definitions  and  Terminology 

This  Section  gives  the  major  definitions  required  for  the  formulation 


of  mathematical  models  in  the  study  of  incentivised  availability  warranties. 

The  models  will  consider  the  payment  by  the  contractor  under  its 
warranty  obligation  as  a  stochastic  process  taking  place  over  the  warranty 
period.  The  models  considered  in  this  report  will  deal  with  situations  in 
which  a  failure  causes  a  payment  "c"  to  be  made  by  the  contractor  at  the 
time  of  failure.  This  may  represent  replacement  or  repair  costs.  The 
time  for  making  repairs  or  replacements  is  not  considered  in  these  models. 

It  can  be  assumed  to  be  relatively  short  with  respect  to  the  full  warranty 
period.  The  present  models  are  intended  to  give  some  basic  insights  into 
the  incentive  warranty  concept  and  it  seems  undesirable  to  include  the  com¬ 
plications  of  repair /replacement  times  in  these  initial  models.  Further 
research  may  well  include  those  times,  providing  additional  insights  into 
the  picture  of  the  availability  profile. 

In  a  stochastic  availability  model  the  number  of  times  an  equipment 
fails  in  the  warranty  period  is  a  discrete  random  variable.  A  general  model 
can  be  made  with  this  point  of  view  and  one  form  of  such  a  model  is  discussed 
in  Chapter  4.  However  for  the  purpose  of  gaining  insight  into  the  incentive 
warranty  concept  such  general  models  are  far  too  complex.  The  mathe¬ 
matical  and  stochastic  features  of  the  general  model  tend  to  overshadow  the 
details  of  particular  interest  in  studying  the  incentive  concept.  Therefore 
the  major  work  of  this  report,  described  in  Chapters  2  and  3  employs  simpli¬ 
fied  models  in  which  at  most  one  or  at  most  two  failures  can  occur  within 
the  warranty  period.  Technically  these  models  may  be  thought  of  as  con¬ 
ditional  models  formulated  under  die  conditional  probabilities  for  the 


number  of  failures.  Since  the  models  are  self  contained  there  is  no  need  to 
include  the  actual  conditioning  features  in  the  details  of  these  simplified 
models.  The  one  or  two  failure  type  models  contain  the  concepts  necessary 
to  study  incentivised  warranties  and  in  many  practical  situations  (though 
certainly  not  all)  at  most  two  failures  may  be  expected  in  a  reasonable 
warranty  period. 

The  case  of  at  most  one  failure  in  the  warranty  period  is  the  relia¬ 
bility  case.  Thus  its  study,  in  Chapter  2,  constitutes  an  extension  of  the 
RIW  to  incentivised  forms. 

Some  of  the  basic  notation  used  in  these  3tudies  will  now  be  given, 
additional  notation  will  be  provided  as  it  is  required. 

The  warranty  period  is  denoted  by  T  .  It  is  a  contract  parameter 
which  specifies  the  duration  of  the  contractors  obligation  for  assuming 
correction  of  failures  under  the  warranty. 

The  warranty  fee  pool,  denoted  by  W(t)  is  an  amount  of  money  that 
the  customer  will  provide  toward  the  repair  replacement  cost  for  a  failure 
at  time  t.  The  fee  pool  forms  employed  in  this  study  are  of  linear  segment 
type.  They  can  most  simply  be  described  in  the  one  failure  (reliability) 
case. 

The  linear  segment  form  in  the  one  failure  case  is: 

Wj(t)  =  aQ  ,  0<t<t' 

=  aQ  +  b(t-t'),  t'<t<T* 


where  a  constant  pool  value  aQ  is  provided  for  part  of  the  warranty  period. 


specified  by  t  ,  and  the  remaining  amount  increases  linearly  at  a  rate 
measured  by  b.  This  type  of  warranty  fee  pool  is  specified  by  three  para¬ 
meters  (contract  parameters):  aQ,  b,  t'  and  represents  a  reasonable 
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form  for  such  a  pool.  It  is  selected  because  of  its  desirability  as  an  actual 
form  rather  than  for  any  mathematical  convenience.  Therefore  it  acts  as 
a  valuable  model  form  in  studying  the  concepts  of  incentivised  warranties. 
Other  warranty  pool  forms  might  be  employed,  for  example  one  based  on 
exponential  type  functions.  However  it  is  felt  that  the  linear  segment  type 
is  best  suited  for  contract  negotiations  and  for  sound  intuitive  understanding 
of  the  features  implicit  in  incentive  warranty  structures.  The  purpose  of 
this  study  is  to  contribute  to  such  understanding  and  this  was  the  reason  for 
selecting  the  linear  segment  form  of  warranty  fee  pool.  Some  calculations 
are  cumbersome  with  this  form  but  it  is  not  clear  that  other  forms  would 
greatly  simplify  calculation.  Of  course  there  is  no  reason  to  limit  the 
warranty  fee  pool  to  linear  segments  and  any  monotonic  non-decreasing 
function  could  be  used. 

For  more  than  one  failure  the  fee  pool  W(t)  becomes  a  stochastic 
process.  Each  time  there  is  a  failure  some  kind  of  adjustment  must  be  made 
in  W(t)  to  account  for  this  occurrence.  The  treatment  employed  in  this  study 
is  as  follows: 

Let  w^(t)  denote  the  warranty  fee  pool  function  after  k  failures.  In 
this  notation  wQ(t)  =  W^t)  defined  above.  If  the  first  failure  occurs  at  tj 
then 

=  0  ,  t<tL 

W^(t)  =  WQ(t^)  -  C  ,  t^  <  t  <  t  #  +  tj 

=  w^tj)  -  c  +  bft-t^t'),  t'+tj<t<T* 

Should  t '  exceed  the  remaining  time  T*-t^  then  w^(t)  has  the  constant  value 
w  <tj)-c  over  the  remaining  period. 
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If  the  second  failure  occurs  at  time  then: 

w2(t)  =  0  , 

-  -  c  , 

=  -  c  +  Mt-^-t'), 


t<t2 

t2  <  t<t'  +  ^ 

/  * 
t  +  t2  <  t  <  T  . 


In  general  a  recurressive  definition  specifies  w^(t)  in  terms  of 
wk  l^k^  w*iere  ^  *be  fbne  of  the  k**1  failure.  In  turn  W(t)  is  defined  in 
terms  of  these  w^(t)  values  and  the  joint  failure  time  distribution  over  the 
warranty  period.  This  general  formulation  will  be  treated  in  Chapter  4. 

The  random  nature  of  the  failure  events  is  determined  by  the  failure 
time  distributions  and  the  number  of  failures.  When  there  is  at  most  one 
failure  the  failure  tune  defines  the  random  characteristics  of  the  model. 
These  studies  make  use  of  two  kinds  of  failure  time  distributions:  uniform 
and  exponential,  for  detailed  analysis.  General  formulations  are  given  in 
terms  of  general  failure  distributions.  The  exponential  distribution  of  the 
form: 

-Xt 


P[T  <  t]  =  1  -  e' 


is  widely  used  in  reliability  work  and  has  been  found  to  represent  many 
actuax  situations.  Its  relative  simplicity  and  wide  applicability  make  it  a  good 
type  of  distribution  to  use  for  some  of  the  present  studies.  In  addition  many 
people  interested  in  reliability  and  related  fields,  such  as  the  one  under  study 
here,  will  be  familiar  with  exponential  type  models. 

The  distribution  used  most  in  the  one  failure  case  (Chapter  2)  is  not 
the  exponential  but  the  uniform  distribution  having  the  simple  form: 

£tw  a  ib  •  r  1 1 5  9 

=  0  ,  elsewhere 

for  the  probability  density  function  for  the  failure  time  T.  Calculations 


using  the  uniform  density  are  rather  complicated,  though  the  exponential  does 
not  greatly  simplify  calculations  as  shown  e.g.  in  Chapter  2.  It  may  be  felt 
that  the  uniform  distribution  does  not  represent  many  actual  situations  of  random 
failure.  However  the  reason  for  using  the  uniform  distribution  as  a  major 
feature  of  the  mathematical  models  for  incentivised  warranties  is  that  their  use 
allows  detailed  study  of  the  major  features  of  incentive  warranties.  Moreover 
actual  failure  situations  can  be  both  approximated  and  bounded  by  the  uniform 
distribution  type  models.  Their  use  directs  attention  to  the  features  of  the 
warranty  sturcture  rather  than  overshadowing  them  in  the  mathematics  of  more 
involved  random  processes. 

The  models  deal  with  the  stochastic  process  A^(t),  the  payment  by  the 
constractor  under  its  warranty  obligations  when  k  failures  occur.  The  process 
is  denoted  by  A(t)  in  the  unconditional  case  where  the  number  of  failures  is 
unspecified,  as  discussed  in  Chapter  4.  Chapter  2  deals  with  A^t),  the  relia¬ 
bility  case  and  Chapter  3  deals  with  A2(t),  the  simplest  case  that  is  of  availa¬ 
bility  type.  Analysis  of  the  models  considers  the  expected  values  E[A^(t)]  and 
variances  Var[A^(t)]  which  are  taken  as  the  quantities  which  represent  the 
significant  features  of  incentive  warranties.  By  considering  the  contract  para¬ 
meters  present  in  the  warranty  fee  pool,  the  parameters  of  the  underlying 
failure  time  distributions,  and  the  form  of  the  process  A^(t)  an  understanding 
of  the  complex  warranty  structure  is  obtained. 

Section  1.  3  Scope  of  Model  Methodology 

The  mathematical  models  employed  in  this  study  are  of  three  general 
types:  at  most  one  failure,  at  most  two  failures,  and  the  general  case  of  un¬ 
specified  number  of  failures.  Each  model  employes  a  specified  form  of  war¬ 
ranty  incentive  fee  pool,  and  probability  distribution  of  time  between  failures. 
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The  study  is  limited  to  a  single  piece  of  equipment  and  its  availability  profile 
over  a  specified  warranty  period.  Correction  of  failure  by  replacement  or 
repair  is  assumed  to  have  fixed  cost  and  to  be  instantaneous.  All  of  these 
conditions  can,  of  course,  be  generalized  and  it  may  be  of  value  to  do  so  in 
subsequent  studies.  In  this  initial  investigation  the  objective  is  not  to  treat 
the  most  general  or  even  the  most  realistic  situations  but  rather  to  study  the 
concept  of  incentivising  warranties  and  consider  their  generalization  from 
reliability  to  availability  form. 

Each  model  provides  an  opportunity  for  analysis  designed  to  show  the 
effect  on  the  payment  process  A^(t)  due  to  postulated  structure  and  magnitude 
of  contract  parameters  defining  the  incentive  fee  pool.  These  effects  must 
necessarily  be  studied  within  postualted  random  behavior  of  equipment  failure. 
The  results  are  intended  to  show  how  incentivising  warranties  motivates  a 
contractor  to  both  assume  such  warranty  obligations  and  provide  a  desirable 
level  of  availability.  The  analyses  of  this  study  constitute  a  background  for 
understanding  the  incentivising  of  availability  warranties.  Though  the  study 
does  not  intend  to  provide  a  detailed  plan  for  structuring  actual  warranties  it 
should  serve  as  an  initial  guide  toward  applications.  The  incentive  fee  pool 
form  used  in  many  of  the  models,  a  linear  segment  type,  seems  reasonable 
for  actual  contract  structuring  and  negotiation.  This  form  is  based  on  a 
structure  that  can  be  intuitively  appreciated  by  contract  negotiators  with  a 
wide  range  of  technical  background.  Moreover  they  are  based  on  a  few  para¬ 
meters  whose  effects  can  be  directly  related  to  both  contractor  and  customer 
goals  through  analyses  similar  to  those  contained  in  this  study. 

An  important  consideration  in  warranties  is  the  value  of  money.  This 
is  because  the  payment  process  A^(t)  will  take  place  in  the  future,  often 
several  years  after  die  original  contract  for  development  and  manufacture  of 
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an  equipment  is  negotiated.  It  has  become  rather  standard  practice  in  studying 
contract  structuring  and  particularly  in  RIW  research,  e.  g.  in  (1),  to  account 
for  the  value  of  money.  The  present  study  is  primarily  intended  as  a  theo¬ 
retical  investigation  of  the  concepts  discussed  above.  The  modifications  to 
true  value  of  money  do  not  effect  these  theoretical  considerations,  though  they 
must  be  considered  when  using  the  model  analyses  as  guides  for  structuring 
actual  incentivised  warranties.  Therefore  most  of  the  present  study  is  car¬ 
ried  out  without  regard  to  value  of  money  considerations.  However  to  provide 
some  indications  of  how  value  of  money  considerations  may  be  included  in  the 
analyses  some  of  the  work  reported  here  does  include  them.  Though  the 
treatment  of  present  value  of  future  expenditure  is  "well  known"  in  certain 
technical  areas  the  basic  ideas,  as  employed  in  the  present  study,  are  included 
here  for  completeness. 

If  an  expenditure  of  amount  s  is  to  be  made  at  time  t  the  present 
value  that  may  be  assigned  to  that  expenditure  is  determined  by  the  discount 
rate  r.  This  rate  is  a  dimensionless  quantity  which  specifies  how  the  value 
of  money  decreases  over  time  due  to  a  combination  of  factors. 

At  time  t  the  expenditure  s  is  made.  Its  value  one  time  unit  (e.  g.  a 

year)  earlier  is  s-rs  or  (l-r)s.  Similarly  if  the  value  k  time  units  earlier  is 

denoted  by  vfc  then  ;  =  (l-r)vfc  so  that  at  the  present  time  (t=0)  the  value  is 

(l-r)ts.  Some  work  dealing  with  present  value  considerations  employes  an 

exponential  form  for  the  discount  rate  (e.  g.  in  (1) ).  In  this  formulation  the 

"rate"  is  most  often  defined  as  a  quantity  y  having  the  dimension  of  inverse 

time.  The  reduction  factor  is  written  as  an  exponential  e”Y  *°  *****  present 

-  Yt 

value  of  an  expenditure  of  amount  s  at  time  t  becomes:  s  e 

More  than  one  future  expenditure  can  be  combined  into  a  present  value 
calculation  as  illustrated  by  the  following  case  of  two  such  expenditures  c1  and 


1-11 


c,  taking  place  at  times  and  t2  respectively.  Note  that  in  the  models 
considered  here  the  discount  rate  is  assumed  to  be  constant.  In  the  kind  of 
applications  under  consideration  in  this  study  the  expenditures  c  ^  and 
may  represent  the  actual  cost  to  the  contractor  under  its  warranty  oblications 
i.  e.  c-Wftj)  and  c-W(t2).  Here  (t^,  t^)  is  a  realisation  of  the  bivariate 
random  variable  (T^,  T2).  The  present  value  of  such  expenditures  (a  random 
variable)  is 

*1  *2 
v  =  c^(l-r)  +  c2(l-r) 

=  cjjl-r)  1  +  (1-r)  -  WftjHl-r)  -  W(t2)(l-r)  c  . 

r  t2~i 

In  this  expression  c|_(l-r)  +  (1-r)  J  gives  the  present  value  of  expenses 

f  1  *2 

when  there  is  no  incentive  fee  pool.  WftjHl-r)  +  W(t2)(l-r)  gives  the 
present  value  of  the  amount  contributed  by  the  fee  pool. 

Though  most  actual  situations  may  be  expected  to  involve  only  a  few 
failures  with  their  associated  costs,  it  is  interesting  to  consider  a  special 
case  that  allows  any  number  of  failures  to  occur.  Such  cases  are  widely 
applied  as  relatively  simple  mathematical  approximations  to  lifetime  avail¬ 
ability  profiles  and  indicate  bounding  values  of  equipment  behavior.  To  pro¬ 
duce  a  useful  model  of  this  kind  it  is  helpful  to  assume  that  the  failures  occur 
at  regular  time  intervals  of  duration  T.  This  is  the  kind  of  model  obtained 
from  so  called  expected  value  techniques  where  T  is  the  mean  time  between 
failures  (e.  g.  (1)  employes  such  an  approach  in  some  of  its  development).  In 
such  a  model  the  present  value  of  all  expenditures  is: 

•  OB 

v  =  c  l  (l-r)kT  -  £  W(kT)(l-r)kT 
k=0  k=0 
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where  it  is  assumed  that  the  first  failure  occurs  at  time  aero.  The  present 
value  of  expenditures  when  there  is  no  incentive  pool,  represented  by  the 
first  infinite  sum  above,  is  equal  to: 


1-(I-DT  ’ 

a  common  form  of  expression  for  the  present  value  of  a  series  of  equal 
expenditures  projected  infinitely  far  into  the  future.  The  second  expression 
is  in  general  complicated  by  the  presence  of  the  stochastic  process  W(t). 
However  it  represents  the  present  value  of  the  incentive  fee  pool  process.  A 
bounding  special  case, which  provides  at  least  some  insight  into  the  full  present 
value  quantity  v  may  be  obtained  by  considering  the  situation  in  which  the 
incentive  pool  is  constant  and  equal  to  a  value  d.  All  this  means  is  that  the 
cost  is  reduced  from  c  to  c-d  at  each  expenditure  and 


v 


c-d 

1  -  (l-r)T 


A  more  interesting  special  case  is  obtained  by  assuming  that  the  fee  pool  values 
increase  to  approach  d.  An  appropriate  form  is  W(kT)  =  d(l-e“kT  where  d  <  c. 


which  is  greater  than  in  the  simple  constant  fee  pool  case,  as  it  should  be, 
representing  the  present  value  of  contractors  expenditure  under  the  warranty 
obligation.  The  exponentially  increasing  fee  pool  provides  a  smaller  contribu¬ 
tion  from  the  customsr. 
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Chapter  2  At  Most  One  Failure  (The  Reliability  Case) 

This  Chapter  defines  and  analyzes  models  representing  incentivised 
reliability  warranties.  At  most  one  failure  is  allowed  in  these  models.  The 
first  model,  given  in  Section  2.  1,  is  very  simple.  It  is  a  deterministic  failure 
time  model,  included  to  provide  a  limit  case  for  the  other  models.  Section  2.  2 
uses  the  linear  segment  type  fee  pool  with  three  cases  of  uniform  failure  time 
distributions:  in  case  1  the  failure  must  occur  during  the  period  of  constant 
pool  value  aQ,  in  case  2  it  may  occur  during  constant  value  a^  or  while  the 
fee  pool  is  increasing  at  constant  rate  b,  the  third  case  deals  with  failure 
occurring  only  within  the  increasing  pool  time  interval.  A  fourth  case  could  be 
considered,  allowing  the  failure  to  occur  toward  the  end  of  the  increasing  pool 
period  or  after  the  warranty  period,  specified  by  time  T  ,  has  ended.  This 
case  does  not  seem  to  contribute  basic  information  of  value  not  already  given 
by  the  other  cases  and  is  therefore  not  included  as  part  of  these  studies. 

The  uniform  distributions  of  Section  2.  2  are  used  to  provide  detailed 
structural  analysis  of  how  the  fee  pool  operates  as  a  motivational  instrument. 
They  are  used  as  tools  of  analysis  rather  than  selected  as  representing  forms 
of  greatest  analytic  convenience  or  representational  accuracy.  The  later  points 
of  view  often  call  for  use  of  exponential  type  failure  distribution  and  these  are 
considered  in  Section  2.  3.  The  respective  utility  of  the  two  types  of  distribu¬ 
tions  is  indicated  within  their  individual  sections  of  analysis.  It  may  be  noted 
that  the  uniform  type  failure  distribution  allows  an  interesting  limit  process 
to  be  formulated,  converging  to  different  values  depending  on  the  point  of  view 
one  assumes.  One  possibility  leading  to  the  deterministic  case  of  Section  2.  1 
as  expected.  However  the  existence  of  different  limit  results  represents  a 
mathematical  formulation  of  certain  qualitative  considerations.  This  is  an 
instance  of  the  trend  in  some  recent  applications  of  mathematics  away  from 
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purely  quantitative  results,  as  for  example  in  the  philosophy  and  application 
of  catestrophy  theory. 

Section  2.4  discusses  the  findings  implicit  in  the  analyses  of  the  other 
sections.  It  also  includes  some  considerations  dealing  with  the  present  value 
of  money.  For  simplicity  the  basic  models  do  not  include  the  concept  of  present 
value  since  they  are  intended  as  theoretical  studies  rather  than  models  of 
actual  warranty  structures. 

All  sections  use  the  reliability  case  fee  pool  of  the  form: 

W(t)  =  aQ  ,  0<  t<  t' 

=  aQ  +  b(t-t/),  t'  <  t  <  T*  . 

The  analyses  of  this  chapter  develop  expressions  for  EtA^t)]  and 
Var[Aj(t)],  the  expected  value  and  variance  of  the  contractor  payment  process 
A^(t)  in  the  case  of  at  most  one  failure.  Expected  values  are  computed  using 
condition  type  logic  of  the  form:  For  a  random  variable  X  defined  over  an 
interval  (O,  R)  with  distribution  function  Fx(x) 

R  R1  R 

E[X]  =  J  x  d  Fx(x)  =  J  *dFx+/  ldFx 
0  0  Rx 

where  0  <  R^  <  R.  Thus  in  dealing  with  the  linear  segment  fee  pool  differing 
payment  forms  are  divided  into  appropriate  segments  for  computing  expected 
values.  The  variance  formulas  make  use  of  the  form: 

VarCx]  =  E[X2]  -  (e[x])2 

properly  divided  to  account  for  segmented  parts  of  the  random  variable 
definition  as  detailed  within  the  appropriate  sections. 

These  basic  ideas  can  be  used  to  formulate  general  formulas  for 
E[A^(t)]  and  VarCA^(t)]  to  be  used  in  the  analyses  of  this  chapter. 
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Here 


AjW  =  0  .  t  <  T 

=  c-W(T),  T  <  t  <  T* 

where  the  random  variable  T  is  the  failure  time  governed  by  a  probability 
density  function  fT(t). 

In  the  following  development  of  basic  formulas  P[  ]  denotes  the 
probability  expression  and  a  vertical  line,  as  usual,  indicates  conditional 
expressions. 

E[Aj(t)]  =  E[A1(t)|t<  T]  PCt<  T ]  +  E[A j(t)  1 1>  T]  P[t>T]  . 

Since  E[Aj(t)  |t  <  T]  =  0  it  follows  that: 

E[Aj(t)]  =  E[c-W(t)|t  >  T]  P[t  >  T],  or  using  the  density  function  fT(t): 

,t  t 

E[Aj(t)]  =  J  [c-W(v)]fT  (vlv  <  t)dv  J  fT(v)dv 
0  0 


I  fT(v) 

However  the  conditional  density  is:  f^,(vl  v  <  t)  =  ■  . . 

J  fT(v)dv 
0 


So  that: 


t 

E[Aj(t)]  =  J  Cc-W(v)]fT(v)dv. 
0 


Considering  the  expression  for  variance  gives: 
VartAjW]  =  E[(Aj(t)  -  ECA^t)])2]  , 
t 

Var[Aj(t)]  =  J  (c-W(v)  -  eCAj^j)  fT(v)dv  , 
0 


Alternatively,  using: 
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VarCAjft)]  =  E^Ajft))2]  -  (eCa^I)])2 
gives  the  formula: 

t 

VarCAjtt)]  =  J  (c-W(v)/fT(v)dv  -  (eCA^)])2  . 

0 

These  formulas  for  expected  value  and  variance  are  employed  in  the  analyses 
of  this  chapter. 

Section  2.  1  Deterministic  Failure 

When  the  failure  occurs  at  a  specified  time  s  with  probability  one  the 
deterministic  failure  model  results.  It  is  the  simplest  case  and  provides  an 
extreme  limit  for  the  other  models,  particularly  the  uniform  case. 

Model  analysis  derives  the  expected  value  and  variance  of  the  contractor 
payment  process  Aj(t).  They  are  expressed  as  functions  of  time  t.  These 
quantities  represent  the  contractors  position  throughout  the  warranty  period 
which  terminates  at  time  T  . 

In  the  deterministic  case  the  derivation  of  ECA^(t)]  and  Var[A^(t)]  is 
immediate,  yielding: 

ECAjft)]  =  0  ,  t<s 

=  c  -  a  ,  t>s  and  s  <  t ' 

O  " 

=  c  -  a  -  b{t-t ,),  T*  >  t  >  s  and  s  >  t ' 

Var[Aj(t)]  =  0  for  all  t. 

The  action  of  the  incentive  fee  pool  is  obvious  in  this  case.  So  long  as 
a<J+b(T*-t/)  <  c  failure  will  cost^the  contractor  money.  It  will  cost  less  when 

it  occurs  further  along  in  the  warranty  period.  This  may  be  considered  the 
classical  situation  from  the  customers  point  of  view.  However  the  cost  re¬ 
duction,  without  an  actual  profit  motive  might  not  be  enough  to  make  such  a  con¬ 
tract  attractive  to  the  contractor.  On  the  other  hand  there  is  no  risk  assumption 


in  this  simple  case.  The  failure  time  is  known  exactly  and  therefore  specifies 
a  known  payment  requirement.  Since  the  opportunity  for  profit  should  imply 
a  degree  of  risk  assumption  on  the  part  of  the  contractor  it  would  be  meaning- 

&  i 

less  to  include  such  a  possibility  (a+b(T  -t  )  >  c)  in  this  simple  case.  When 
the  failure  time  is  uncertain,  as  it  most  often  is  in  practice,  more  possibili¬ 
ties  exist. 

Though  this  simple  case  is  meant  only  as  a  limit  consideration  and  an 
introduction  to  the  analysis  approach,  it  can  be  thought  of  as  a  practical  model 
in  one  sense.  When  preventive  replacement  is  carried  out  at  assigned  times 
that  are  short  enough  to  essentially  insure  no  failures,  the  model  represents 
the  situation  by  interpreting  deterministic  failure  as  scheduled  preventive 
replacement. 

Section  2.  2  Uniform  Failure 

In  this  section  the  failure  time  is  determined  by  a  uniform  distribution 
having  density  function: 

V*)=,-Tr'  r-‘^s 

=  0  ,  elsewhere. 

The  analysis  differs  in  details  of  calculation  and  results  as  the  uniform  density 
is  placed  so  as  to  completely  proceed,  include,  or  completely  follow  the  incen¬ 
tive  fee  parameter  t' ,  Each  of  these  three  cases  will  be  treated  individually. 

The  calculations  use  the  general  formulas  given  at  the  beginning  of  this 
chapter. 

Case  1  s  < 

Direct  calculation  yields: 


E[Aj(t)]  =  0  ,  t  <  r 

t 

=  J  <c-ao>ibdv  = 

.  r 

s 

=  J  (c-ao)ndv=(c-ao,(  T^>t>3 

r 

E[Aj(t)]  =  0  ,  t  <  r 

=  (c-ao)2^TF'  '^l3 

=  (c-aQ)2  ,  T*  >  t  >  s 

Var[A.(t)]  =  (c-a  )2  i^Sli3^).  f  r<t<  s 

(s-r) 

=  0  ,  elsewhere. 

An  alternative  form  of  these  results,  particularly  useful  for  limit  studies,  is 
obtained  by  the  change  of  variables:  r  =  s-e,  t  *  s  -  P  e,  0  <  P  <  1.  Only 
the  range  r  <  t  <  s  need  be  considered  (the  other  values  are  properly  assumed 
at  the  end  points  of  this  range). 

ECAj(t)]  =  (c-ao)(l-P) 

VarCAjW]  =  (c-aQ)2  P(l-P). 

As  the  width  of  the  uniform  distribution  changes  the  expected  value  and  variance 
of  the  payment  process  A^(t)  are  not  effected.  They  depend  only  on  where  t 
is  located  within  the  range  of  definition  of  the  uniform  distribution.  As  t  tends 
to  the  end  value  of  the  distribution,  s,  the  parameter  {3  becomes  zero  giving 
the  limit  values  obtained  for  the  deterministic  occurance  of  a  failure  at  time 
t=s  in  Section  2.  1.  These  qualitative  limit  considerations  will  be  discussed 
further  under  case  3. 

2.6 
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Case  2  r  <  t  <  s 

It  is  only  necessary  to  consider  the  mean  and  variance  expressions  for 
r  <  t  <  s  and  obtain  the  proper  expressions  for  other  values  of  t  from  the  end 
point  values  of  this  range.  Thus  the  following  calculations  asstime  the  range 
r  <  t  <  s. 

t 

ECAjM]  =  f  (c-ao)^dv  =  <c-ao)££  • 

r 

y  t 

ECAj{t)]=j  (c.ao)^+f  (c-ao-b(v-t/))^X  ,  t'<t<s 


,  x  t-r  b  (t-t/)^ 

'  o'  s-r  2  (s-r) 


,  t  <  t  <  s 


E[Al(t|].  fc-ao)-|%4  .  t>,. 


(s-r) 

The  last  expression  is  obtained  by  letting  t=s  in  the  previous  expression.  In 
the  remaining  calculations  it  is  convenient  to  introduce  the  quantity 
a  =  c_a0  +  bt*  as  a  combination  of  fee  pool  parameters  and  cost  that  occurs 
throughout  much  of  the  work.  Thus  the  complete  result  is: 

E[Aj(t)]  =  0  ,  t<  r 

=  (c-ao}  7T§  » 

a  (c-ao)  -  \  ,  t>*. 

Calculations  also  yield: 
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t  <  r 


E[Aj(t)]  *  0 


/  \2  t~r 

=  (c-a  )  - 

'  o'  s-r 


r  ^  t  <  t 


2„  '3  3. 


•  c=-s)2irf  - ^  («2  -  »■»(*** '))  - b  %:-1y 1  ■  ‘'i1-3 
3  <c-V2 Tif  - (»2  - « -  bW  •  'i* 


These  values  lead  to  the  following  expressions  for  the  variance: 


VarCA^t)]  =  0 


,2  (t-r)(s-t) 
=  (c-a  )  5 —  M 

°  (s-rr 


,  t<  r 


»  r  ^  *  £  t 


For  t  >  s  the  expression  is  obtained  from  the  last  formula  with  t-s. 

Case  3  r  >  t' 

In  this  case  there  is  no  probability  of  failure  during  the  constant  fee 
pool  period  t  <  t'  so  that  there  is  an  increasing  incentive  throughout  the  entire 
range  of  definition  of  the  uniform  failure  time  distribution.  This  is  the  most 
interesting  of  the  three  cases  considered,  and  a  more  extensive  analysis, 
including  some  limit  studies  will  be  given  for  it.  The  basic  results  for  this 
case  follow  from  application  of  the  general  formulas,  yielding  the  following: 


ECA^t)]  =  0 


.  t  <  r 


(•-TO  * 


-a.Sfef£L 


,  s  <  t  <  T 


where  a  =  c-a  +  bt  . 

o 
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I 


2 

E[Aj(t)]  *  0  ,  t  <  r 

2 

=  [a2  -  ba(t+r)  +  (t2  +  tr  +  r2)  |  ,  r  <  t  <  s 

=  a2  -  ba(s+r)  +  (s2  +  sr  +  r2)  ,  3  <  t  <  T* 


These  values  yield  the  expression  for  variance: 
V-.rCAj(t)]  =  0  ,  t<t 


■•t~£-yjra2  -  ab(ttr)]  (s-t)  +  jy[-3t3  +  (4s  -  7r)t2 
(s-r)  |L  . 

+  (4sr  -  r2)(t+r)j|  ,  r  <  t<s 


•yj  (s-r)2  »  3  T 


The  end  range  values,  s  <  t  <  T  can  be  obtained  in  an  alternative,  simple  way 
using  the  properties  of  the  uniformly  distributed  failure  time  T.  This  simple, 
direct  calculation  is  included  here  for  its  interest  and  to  illustrate  the  underly¬ 
ing  role  of  the  failure  time  process. 

A^t)  =0  ,  T>t 

Ax(t)  =  a  -  bT  ,  T  <  t 

E[Aj(t)]  =  E[a]  -bE[T]  for  s  <  t  <  T* 


the  relation  between  the  running  time  value  t  and  the  random  time  to  failure 
T  are  such  as  to  cause  this  relation  to  only  apply  in  the  specified  time  range. 
Since,  for  the  uniform  random  variable  T,  E[T]  =  —j^,  and 

Var[T]  =  ■ ,  the  following  results  are  obtained 


ECA^t)]  =  a 

VarCAjCt)]  =  VarCa]  +  b2  Var[T] 


as  given  above  from  the  direct  calculation  approach. 

An  alternative  form,  useful  for  limit  considerations  is  obtained  by 
making  the  change  of  variables:  r=s-e  and  t=s-0  e  where  0  <  0  <  1.  Only  the 
mid-range  expression  for  r  <  t  <  s  will  be  presented  since  the  other  values 
may  be  obtained  from  these  as  special  cases. 

ECA^t)]  =  [a  -  bs  +  b(l-0)  f  j  (l-P) 

Var[Aj(t)]  =  (l-0)|[a2  -  ab(2s  -  (l+p)e)]p 

2  -i 

+  [(303  +  702  +  P+l)  e2  -  12  s  0  (P+1)  e  +  12  0  s2_i 

In  this  formulation  e  represents  the  interval  in  which  a  failure  takes 
place.  It  may  be  considered  very  small  in  order  to  yield  an  interesting  special 
form  that  allows  very  detailed  consideration  of  how  the  payment  relates  to  the 
time  at  which  it  is  considered.  The  simplest  expression  corresponding  to 
small  e  is  to  set  e=o  in  the  above  formulas  resulting  in: 

EL[Ax(t)]  =  (a  -  bs)  (1  -  0) 

Var^A^t)]  =  (1  -  0)  0(a  -  bs)2 

where  the  subscript  L  denotes  the  extreme  value  for  e=0. 

It  is  seen  that  E[Aj(t)]  and  Var[Aj(t)]  do  not  tend  to  a  limit  in  the 
mathematical  sense  as  e-0.  The  limit  type  expressions  give  different  values 
for  various  0.  From  a  limit  povnt  of  view  E[A^(t)]  and  Var[A^(t)]  depend  on 
two  variables:  the  width  e  and  the  position  of  t  within  the  distribution,  de¬ 
termined  by  0.  In  a  gross  sense,  e  -0  forces  t-s  for  all  values  of  0.  How¬ 
ever  the  above  results  show  that  the  limit  depends  on  the  values  assumed  for 
0.  If  0=0  the  evaluation  time  t  is  held  at  s,  corresponding  to  the  values 
obtained  for  the  deterministic  situation  of  Section  2.  1.  In  this  case  one  is 
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considering  the  extreme  form  of  the  payment  function  after  the  failure  has 
occurred,  there  being  no  density  of  probability  left  over,  after  time  t,  to  allow 
the  failure  to  occur  later.  On  the  other  extreme  with  P  =  1  the  time  is  considered 
before  the  probability  density  starts  so  that  the  failure  could  not  have  occurred 
and  of  course  E[A^(t)]  =  0  in  this  case.  In  both  of  these  extreme  j3  cases 
the  variance  is  zero,  as  it  should  be  because  there  is  no  variation  of  the  failure 
event.  In  the  P=0  extreme  the  event  is  sure,  in  the  jS  =  l  extreme  the  event  does 
not  occur.  Between  these  extremes  the  evaluation  time  is  within  the  probability 
density  (strictly  speaking  e  is  small,  having  little  effect,  but  not  zero;  however 
the  simpler  formulas,  with  e=0  are  useful  as  approximations  to  the  more 
general  expressions).  In  such  cases  there  is  a  chance  that  the  failure  has 
occurred  by  time  t.  There  is  also  a  chance  that  failure  will  occur  after  time  t. 
The  expected  payment  increases  with  t  because  the  chance  of  failure,  the  event 
requiring  a  payment,  increases  with  t  corresponding  to  p-0.  On  the  other 
hand  the  variance  is  greatest  in  the  mid-range,  falling  off  at  each  extreme 
because  the  actual  occurrence  of  the  failure  event  is  most  uncertain  in  the  mid¬ 
range. 

These  are  qualitative  considerations  expressed  in  mathematical  terms 
and  measurable  to  a  degree  by  the  quantitative  parameter  p.  In  applications  they 
may  give  some  guidelines  on  the  level  of  resources  (represented  here  by  pay¬ 
ments)  that  should  be  made  available  at  a  time  t  within  a  relatively  short  period 
wherein  its  expenditure  will  be  required  (for  example:  such  situations  can 
arise  from  scheduled  replacement/repair,  or  high  stress  use  of  equipment  at  a 
specified  time). 

* 

A  particular  instance  of  these  extreme  forms  takes  place  when  s=T  , 
the  end  of  the  warranty  period.  This  may  be  called  the  One-Hoss  Shay  (OKS) 
case  after  "The  Deacon's  Masterpiece,  or,  the  Wonderful  'One-Hoss  Shay'", 
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by  Oliver  Wendell  Holmes,  in  which  the  Shay  operates  perfectly  until,  all 
at  once  it  completely  fails.  The  formulas  for  this  situation,  an  ideal  from 
the  point  of  view  of  the  incentives ed  warranty,  are: 

EOHS[Al(t)]  =  («-bT*)(l-f3) 

VarOHS[Al(t)3  =  te-bT*)2- 

For  completeness,  and  as  a  check  on  some  of  the  calculations,  it  is 
interesting  to  compare  results  for  case  2  and  case  3.  When  t^r  both  cases 
should  yield  the  same  results.  This  is  the  case  with: 

E[Aj(t)]  =  [c-aQ  -  ^(t-r)  J  where  t'=r 

in  the  fee  pool  form. 

2 

The  common  value  of  E[A^(t)]  is  -  ab(t+r)  +  ^-(t^  +  tr  + 

and  since  the  expected  values  and  expected  values  of  the  squares  are  equal  it 
follows  that  the  variance  expressions  are  also  equal  as  desired. 

It  may  be  noted  that  an  interesting  situation  is  obtained  by  setting  t^O. 
This  may  be  called  the  increasing  incentive  type  of  warranty  having  no  initial 
period  of  constant  (possibly  zero)  fee  pool.  Though  an  important  special  case 
for  practical  consideration  it  presents  no  special  form  for  the  expressions 
ECAjft)]  or  VarCAjft)]  other  than  to  observe  that  in  this  case  a=c-aQ. 

Section  2.  3  Exponential  Failure 

In  the  models  described  in  this  section  the  incentive  fee  pool  is  the 
same,  linear  segment,  type  used  in  previous  sections.  The  failure  time  proba 
bility  density  function  has  the  form: 

fT(t)  =  Xe‘U  ,  t>0 

=  o  elsewhere  . 
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This  kind  of  failure  time  distribution  gives  a  very  different  kind  of  model 
than  those  resulting  from  uniform  failure  time  distributions.  The  structure  is 
simpler,  requiring  less  cases  for  consideration  and  may  reflect  some  failure 
situations  more  realistically  than  the  uniform  densities  could.  However  the 
exponential  failure  form  does  not  allow  a  detailed  study  of  the  payment  process 
such  as  can  be  obtained  from  placing  uniform  densities  of  selected  width  at 
various  positions  in  the  warranty  period.  The  two  kinds  of  distributions  com¬ 
plement  each  other  in  the  study  of  the  incentive  warranty  concept.  Though 
other  types  of  failure  distributions  might  be  employed  for  practical  application 
only  these  two  are  given  detailed  analysis  in  the  report.  Some  discussion  of 
general  formulations  is  given  in  Chapter  4. 

One  difference  from  the  uniform  cases  of  Section  2.2  is  that  there  is 
always  some  non-zero  probability  of  the  failure  event  occurring  after  the  end 
of  the  warranty  period,  designated  by  T*.  Therefore  the  expected  values  are 

computed  as  functions  of  the  evaluation  time  t,  never  yielding  the  full  expecta- 

£ 

tion  over  all  non-zero  failure  events.  The  payment  process  A^t)  ends  at  T 

$ 

and  by  setting  it  to  zero  for  t>  T  one  can  interpret  the  expected  value  formulas 

)jc 

(including  variance)  given  below  as  being  completed  by  the  time  t  reaches  T  . 

In  this  interpretation,  though  the  expected  values  are  completed  the  probability 

of  zero  payment  is  not  zero  but  is  equal  to  the  probability  density  remaining 

after  the  expiration  of  the  warranty.  Denote  the  probability  of  zero  payment  by 

P„  then: 
o 

Po  =  J  Xe-Xtdy  =  e-*T* 

T 

this  is  a  useful  quantity  in  the  exponential  failure  models.  A  similar  quantity 
can  be  defined  in  any  model  assigning  non-zero  probability  density  to  time  events 
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after  T  .  In  studying  uniform  distributions  this  situation  was  not  considered 
in  the  present  study  since  it  is  reasonably  represented  in  the  present  model. 


ECAjft)]  =  J  (c-aQ)  Xe‘Xvdv  ,  0<t<t‘ 


=  J  (c-aQ)  X  e"Xvdv  +  f  (a-bv)  X  e  "^vdv,  t'  <t<T^ 


These  values  become: 

ECAjft)]  =  (c-aQ)  (l-e'U)  ,  0<t<t' 

=  (c-aQ)  (1-e  )  +  b(t-t  )e  +  ■£■  (®  "e  ),  t<t<T 

It  may  be  observed  that  the  expressions  for  E[A^(t)]  agree  at  t-t' ,  as  they 
should. 

The  term  (c-ao)(l-e"Xt)  is  the  expected  contractor  cost  with  only  a 
constant  incentive  aQ  provided  by  the  customer.  The  remaining  expressions 
represent  the  contribution  of  the  increasing  incentive  pool  values  to  the  con¬ 
tractors  payment  process  under  its  warranty  obligations.  If  the  incentive  is  to 
be  effective  these  terms  should  be  negative  for  all  t>t*  thereby  reducing  pay¬ 
ment  that  may  occur  at  time  before  the  evaluation  time  t.  The  sign  of  these 
terms  is  determined  by  the  sign  of  the  expression: 

X(t-t')  +  l-e'X(t 

which  may  be  shown  to  be  negative  for  all  t  as  required.  Moreover  the 
expression  becomes  increasingly  negative  as  t  increases  as  it  should  for 
proper  incentive  sturcture. 

2 

The  expression  for  E[A^(t)]  is  as  follows: 
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E[A2<t)]  =  (c-aQ)2  (l-e'Xt)  ,  0<t<t' 

=  (c-ao)2  <l-e_Xt  )  -  [(a-bt-  x)2  +  ^7]  e‘Xt 

X 

This  expression  indicates  a  complicated  form  for  Var[Aj(t)]  in  the  exponential 
case.  Because  of  this  complexity  the  variance  is  not  formulated  in  this  study. 
However  the  special  case  where  t  =0,  with  no  constant  fee  pool  period  may  be 
considered.  This  is  the  true  incentive  case  of  most  interest  and  the  results, 
while  not  particularly  simple,  are  reasonable. 

In  the  case  t/=0,  a=c-aQ  and  there  is  only  one  range  of  definition, 

$ 

(0,  T  ),  for  the  functions  developed. 

ECA^t)]  =  <c-aQ  -  £)  -  (c-aQ  -  £  -  bt)e_Xt 

E[Aj(t)]  =  (c-aQ  -  x)  +  -  [(C"a0  -  X  "  bt)2  +  “z]  e'U 


.  b  ...2  -Z\t 

-  (c-aQ  -  ^  -  bt)  e 

Consider  the  end  value,  at  t=T*,  which  measures  the  expected  payment 
if  a  failure  occurs  within  the  warranty  period.  For  simplicity  the  case  t  =0 
is  presented: 

$  * 

E**E[A1(t)]  =  (c-aQ  -  x)U-«’XT  )  +bT*  e'XT 


If  instead  of  a  warranty  period  of  duration  T  a  lifetime  maintenance  was 
assumed  by  the  contractor,  with  the  same  kind  of  incentive  pool  and  one  failure 


2. 15 


assumed,  the  resulting  expected  payment  will  be: 


®  b 

E  =c-aoT 


Of  course  it  is  unrealistic  to  have  a  continually  increasing,  unbounded  fee 
pool.  However  such  an  assumption  is  useful  here  for  simplicity  and  is 
appropriate  because  the  probability  of  failure  decreases  strongly  as  time 
increases.  The  combined  effect  is  for  a  reasonable  expected  value,  as  given 
above. 

The  quantity  E°°  -  E*  represents  the  difference  in  payments  by  the 
contractor  under  lifetime  maintenance  and  warranty  type  contracts.  Of  course 
the  major  assumptions  are  at  most  one  failure,  and  the  same  incentive  fee 
pool.  Though  somewhat  unrealistic  as  practical  contract  assumptions,  this 
illustration  gives  an  idea  of  one  methodological  approach  for  comparing  dif¬ 
ferent  maintenance  concepts. 

e"  -  E*  =  Cb  T*  -  (c-aQ  -  £)]  e'^T  , 

and  depending  on  the  parameter  values  different  situations  may  be  considered. 
A  case  of  some  interest  occurs  when  the  factors  effecting  the  payment:  failure 
distribution,  and  fee  pool  form,  are  such  as  to  give  the  same  expected  pay- 

CO  ♦ 

ment  with  both  methods .  This  occurs  when  E  -  E  =  0  or: 

,  *  b 

bT 

A  result  such  as  this  indicates  the  interrelation  of  contract  parameters, 
aQ,  b,  T  ,  replacement  cost  c,  and  failure  distribution  parameter  To 
the  extent  that  such  expressions  could  be  developed  (even  in  tabular  form  for 
complicated  situations)  to  represent  actual  situations,  they  would  be  useful 
negotiation  tools. 


Section  2.4  Discussion,  Including  Present  Value  of  Mone 
The  payment  by  the  contractor  under  an  incentive  warranty  is  a 
stochastic  process  of  simple  type,  defined  as: 


=  0 

T  >  t 

*  c-ao 

T  <  t'. 

T  <  t  <  T 

s  c  -aQ  +  b  t  -  bT  , 

T  >  t\ 

T<  t  <  T 

=  0 

T  >  T  . 

The  actual  payment  is,  in  any  case,  a  numerical  value  determined  by  the  time 
at  which  the  random  failure  event  takes  place,  denoted  by  the  random  variable 
T.  Thus  A^(t)  is  a  step  function  of  t  with  the  step  location  determined  by  T 
and  the  step  height  determined  by  the  fee  pool  value  at  the  time  the  step  occurs. 

In  motivating  the  contractor  to  increase  the  time  to  failure  by  providing 

an  increasing  fee  pool  there  is  a  danger  that  the  contractor  will  desire  a  failure 

toward  the  end  of  the  warranty  period  so  as  to  receive  incentive  fee.  This 

undesirable  situation  can  be  prevented  by  selecting  contract  parameters  so  that 

c-a  +  bt'  -  bT*  >  0  so  that  any  failure  will  actually  cost  the  contractor  money, 
o 

There  is  another  point  of  view  however  that  considers  the  provision  of 
an  incentive  warranty  fee  pool  insufficient  inducement  for  the  contractor  to 
agree  to  a  warranty.  It  is  true  that  its  risk  of  payment  is  reduced,  but  with 
the  condition  of  possible  negative  payment  the  warranty  still  represents  a 
(randomly  specified)  cost  to  the  contractor.  Therefore  there  is  a  problem 
with  the  basic  approach:  in  order  to  avoid  the  counter  motivation  of  wanting  a 
failure  the  chance  for  obtaining  profit  by  agreeing  to  a  warranty  is  lost.  One 
way  to  combine  the  desired  motivation  is  to  agree  to  an  award  fee  of  some  kind 
if  no  failure  occurs  by  the  end  of  the  warranty  period.  A  reasonable  value  for 
such  an  award  fee  might  be  aQ  -  b  t '  +  b  T*,  the  value  reached  by  the  incentive 


fee  pool  by  the  end  of  the  warranty  period.  If  availability  is  a  major  concern 
to  the  customer,  with  savings  in  money  a  secondary  consideration,  then  a 
substantial  award  fee  suggests  itself.  In  the  combined  incentive  warranty  pool, 
award  fee  concept  the  contractor  is  motivated  to  assume  the  warranty  obliga¬ 
tion  by  the  chance  of  profit  and  motivated  to  design  for  long  operating  times 
(in  the  at  most  one  failure  case). 

In  the  previous  sections  of  this  chapter  the  payment  process  A^(t)  is 
characterized  by  the  expected  payment  up  to  time  t  and  by  the  variance  of  the 
process.  These  expressions  may  be  useful  for  planning  before  the  start  of 
the  warranty  period  and  as  management  decision  tools  during  the  warranty 
period.  Of  course  once  the  failure  occurs  the  payment  is  known  exactly  and  the 
process  Aj(t)  ceases  to  represent  the  payment.  The  effectiveness  of  the 
models  and  their  formulas  for  expected  value  and  variance  for  planning  and 
control  purposes  depends  to  some  extent  on  the  time  scales  involved.  When 
the  warranty  period  is  long  the  change  in  expected  payment  over  time  can  be  a 
valuable  planning  tool,  particularly  when  the  necessary  payments  under  the 
warranty  are  large. 

The  research  described  in  this  report  is  directed  toward  developing  a 
methodology  for  studying  the  concepts  of  incentivising  warranties.  However 
some  indication  of  the  application  of  this  methodology  can  be  found  in  the  fol¬ 
lowing  discussion.  As  an  example  case  the  third  case  of  Section  2.  2  is  con¬ 
sidered.  In  this  case  the  uniform  distribution  is  used  with  r  >  t' ,  so  actual 
payment  must  occur  in  the  purely  increasing  incentive  region  of  the  fee  pool 

The  variance  is  a  cubic  function  of  t  with  complicated  coefficients. 

It  is  best  studied  by  numerical  calculations  (graphic  plots)  which  are  not 
included  as  part  of  the  present  theoretical  research.  However  the  expected 
value  function,  denoted  by  E(t)  *  ELA^t)]  below,  is  a  relatively  simple 
quadratic  function  of  t  of  the  form: 
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where 


This  function  is  a  parabola  in  the  ^t,  E(t)J  plane  with  vertex  at  (C,A),  cros¬ 
sing  the  time  axis  at  r  and  -  r.  As  discussed  above,  a  reasonable  condi- 

tion  for  the  incentive  warranty  parameters  is  that  a  -  b  t  >  0  for  all  t  <  T 
a  * 

which  implies  that  ^  >  T  in  any  warranty  situation  satisfying  this  condition. 
Under  this  assumption  the  function  E(t)  is  increasing  over  the  portion  of  the 


warranty  period  in  which  the  failure  takes  place,  the  interval  (r,  s)  with 
♦ 

s  <  T  .  The  maximum  value  is  reached  when  t=s  and  from  then  on 

E*(t)  =  a  -  for  s  <  t  <  T*.  This  value  is  easily  shown  to  be  less  than 

the  maximum  of  the  parabolic  function  E(t)  continued  beyond  t=s.  Therefore 


the  relations  between  the  various  time  events,  and  expected  values  are  as 

$ 

shown  in  Figure  2.  1,  where  E(t)  is  the  quadratic  function  and  E  (t)  the  actual 

$ 

expected  value  function,  E  (t)  =  E(t),  r  <  t  <  s. 


Figure  2.  1  Expected  value  function  for  Case  3  Section  2.  2. 
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The  actual  payment,  as  distinct  from  its  expected  value  function,  is 

less  when  the  time  of  the  failure  event  is  great.  The  least  payment  is  at 
*  * 

t  =  T  in  amount  a  -  b  T  and  the  greatest  payment  is  at  t=0  in  amount  a. 
However  in  the  case  under  consideration  here  the  failure  must  occur  in  the 
interval  (r,s)  with  t  <r<s<T  so  the  maximum  payment  in  this  model  is 
a  -  br  and  the  minimum  payment  is  a  -  bs.  The  average  of  these  values  equals 
the  value  of  the  expected  value  function  for  s  <  t  <  T  .  Thus  the  contractor 
may  pay  as  little  as  a  -  bs  or  as  much  as  a  -  br  in  a  particular  situation.  If 
the  same  kind  of  contract  was  repeated  many  times  the  average  payment  over 
all  would  be  1  -  ^  (s+r)  per  contract.  This  average  value  is  the  value  actually 
paid  when  failure  occurs  at  the  average  time  (r+s)/2.  Actual  payment  is  con¬ 
trasted  with  the  expected  value  function  in  Figure  2.  2  for  the  range  (r,s)  in 
which  the  actual  payment  must  be  determined. 


In  Figure  2.2,  Y(t)  represents  the  actual  payment  if  the  failure  event  occurs 
at  time  t,  i.  e.  T=t.  Though  these  are  events  of  probability  sero  due  to  the 
continuous  distribution  assumed  for  the  random  variable  T,  the  curve  Y(t) 
still  serves  as  an  indication  of  what  the  costs  will  be  should  failure  occur  at 
time  t. 
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An  alternative  consideration  can  be  made  utilizing  cumulative  probabil¬ 
ity  statements  about  the  amount  of  payment.  For  example  if  the  random 
variable  Y  is  used  to  denote  payment  then  P[T  <  t]  =  P[Y  <  S  -  bT].  The 
probability  distribution  yields: 

P[Y<a-b  T3=ifi. 

Therefore  if  it  is  desirable  to  keep  the  warranty  payment  below  a  value  a  -  bT 
the  equipment  must  be  designed  to  have  a  time  to  failure  exceeding  T.  The 
strength  of  this  desire  may  be  measured  by  the  probability  •  In  this  situa¬ 
tion  the  contractor  can  spend  money  on  improving  the  reliability  so  that  failure 

occurs  after  time  T  or  he  can  take  the  risk,  measured  by  1 - - - ,  that 

1  s-r  s-r  ’ 

the  failure  will  occur  after  T  with  the  (r,s)  uniform  failure  design.  These 
considerations  illustrate  the  interplay  between  the  incentive  warranty  methodology 
and  equipment  reliability  in  which  various  management  decisions  and  design 
tradeoffs  exist. 

In  all  of  the  discussion  above  the  present  value  of  future  expenditures 

of  money  was  not  considered.  The  purpose  of  the  research  in  this  report  is  to 

introduce  the  concepts  and  model  methodology  for  incentive  warranties.  In 

such  considerations  it  is  not  necessary  to  include  the  actual  value  of  money  (or 

resources)  such  details  falling  into  the  application  area  based  upon  the  theoretical 

developments  of  this  study.  However  to  illustrate  how  the  value  of  money  can  be 

introduced  into  the  methodology  an  example  will  now  be  considered. 

T 

Using  discount  rate  d  a  discount  factor  of  (1-d)  is  applied  to  an  ex¬ 
penditure  of  money  (resources)  at  time  t  to  yield  the  present  values  of  that 
expenditure,  as  discussed  in  Section  1. 3.  This  consideration  can  be  applied  to 
the  payment  itself  or  to  its  expected  function.  Alternatively  various  "expected 
value"  approaches  can  be  taken  to  give  (simpler)  approximate  results.  Such 


expected  value  techniques  are  widely  used  in  the  analysis  of  random  models. 
Therefore  it  is  interesting  to  consider  several  formulations  of  the  present 
value  analysis.  The  uniform  distribution*  Case  3  with  t7  <  r  will  continue  to 
serve  as  the  example  model. 

The  present  value,  considered  as  a  random  available,  has  the  form 
T 

V  =  (a-bT)(l  -  d)  ,  where  T  is  the  uniformly  distributed  random  variable 

in  the  interval  (r,s).  This  formulation  shows  the  disadvantage  of  a  discount 

factor  of  the  given  type  for  calculations  such  as  expected  value.  Therefore 

•lit  -h 

an  alternative  discount  factor  e  is  used,  where  e~  =  1  -  d  relates  the  two 
discount  rate  forms. 

Now  the  present  value  is  the  random  variable: 

V  =  (a  -  bT)  e'hT. 

Direct  calculation  yields: 

E[v]  =  ssb y  [(bs  Hr*/  «'hs  +  (a- bt  -  S  *'hr] 

This  is  to  be  compared  with  the  expected  payment  without  discounting 
E[Y]  =  a  -  Y(s+r)*  '$y  expanding  the  exponential  factors  in  E[v]  one  can 
write: 

-  IT7T [(*4*  +^)(l-h*  +  !Tr  * °<b3>) 


+  -  ^f)(‘  -hr  +  ^T-  +  0(h3>)] 

In  this  form  it  is  possible  to  obtain  the  limit  as  the  discount  rate  goes  to  sero. 

The  result  is  as  expected: 

lizn  E[V]  =  E[Y]  . 
h-0 


At  the  other  extreme  for  total  discounting  the  discount  rate  d=l  with 


r 
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corresponding  condition  obtained  by  taking  the  limit  as  h  becomes  unbounded. 

This  clearly  yields  lim  E[v]  =  0  as  expected. 
h-*» 


For  completness  two  expected  value  type  models  will  be  given  next. 
They  illustrate  techniques  that  are  sometimes  used  and  may  give  results 
ranging  from  correct  to  approximate  to  totally  wrong  and  misleading. 
Expected  value  model  A: 


Ea  =  E[a-bT]  e‘hE^T^ 


ea=  L« -7  <•+*>]  e 


■\  <.+r> 


Expected  value  model  B: 

Eb  =  E[a-bT]  E[e'hT] 


*  (I <8+c>  - «;  *-h’  +  (« -  T  <8«>)  «‘hr 

It  can  be  seen  that  neither  of  these  procedures  give  the  value  E[v]  and  the 
question  of  their  value  as  approximations  would  require  numerical  comparison 
of  the  forms  as  functions  of  h.  However  it  is  interesting  to  note  that  both 
expected  value  models  give  the  non -discounted  expected  value  in  the  limit  as 
h“*0  and  both  tend  to  zero  as  h  becomes  unbounded.  Thus  their  gross  be¬ 
havior  follows  the  correct  value  E[v]  in  the  extreme  cases.  Further  investi¬ 
gation  of  the  value  of  these  quantities  as  useful  approximations  to  ElV]  is 
outside  the  scope  of  the  work  covered  by  this  report.  It  should  be  observed 
that  one  would  expect  that  E[V]  <  eCy]  for  all  values  of  h,  equality  holding 
only  when  h=0.  For  the  expected  value  models  to  be  at  all  useful  they  too 
should  satisfy  this  condition  with  respect  to  the  undiscounted  payment  expres¬ 
sion.  Analytic  varification  of  these  conditions  seems  excessively  complicated 
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and  numerical  checks  would  be  the  most  direct  way  of  varifying  (or  disproving) 
the  relation.  Such  numerical  work  is  outside  the  scope  of  this  initial,  largely 
theoretical,  conceptual  study. 


Chapter  3  At  Most  Two  Failures 

Section  3.  1  Formulation  of  the  Model 

When  more  than  one  failure  can  occur  many  new  features  may  be 
introduced  into  the  warranty  payment  model.  The  material  in  this  report 
represents  an  initial  formulation  of  the  incentive  warranty  structure  and  as 
such  it  treats  rather  simple  models.  The  cost  of  repair  is  considered  a  con¬ 
stant  c  for  each  failure,  no  consideration  is  given  to  time  to  repair /replace, 
nor  are  logistic  support  features  such  as  preventative  maintenance  included. 
Moreover  the  forms  of  the  incentive  pool  and  failure  time  probability  are 
simple.  This  chapter  considers  the  situation  in  which  at  most  two  failures 
can  occur.  The  model  is  considered  from  various  points  of  view  and  the  pay¬ 
ment  process  A^(t)  is  taken  as  the  mathematical  model  for  the  incentive  war¬ 
ranty  structure.  The  expected  value  and  variance  of  A£(t)  are  calculated  as 
being  the  representational  quantities  of  interest.  However  it  should  be  noted 
that  one  could  also  consider  the  probability  that  A^Jt)  has  some  particular 
property  e.  g.  P[A2(t)  <  k]  might  prove  to  be  of  value  in  the  study  of  process 
A2<t).  The  present  report  does  not  work  with  such  probability  values,  re¬ 
stricting  itself  to  development  of  the  representational  moment  expressions. 

In  this  section  the  basic  assumptions  will  be  given,  the  mathematical 
model  formulated,  and  general  expressions  for  expected  value  and  variance 
stated. 

To  model  the  incentive  warranty  process  the  most  basic  quantity 
required  is  the  form  of  the  incentive  fee  pool.  In  the  case  of  more  than  one 
failure  there  is  a  complication  with  the  idea  of  incentive  payment.  If  the 
maximum  payment  is  less  than  c,  the  cost  of  a  single  repair  then  there  is  no 
difficulty.  However  such  a  value  for  incentive  may  not  be  considered  sufficient 


motivation  by  the  contractor  for  assuming  the  warranty  obligation.  Therefore, 
in  the  at  most  two  failure  case  it  seems  reasonable  to  have  the  maximum 
exceed  c  but  be  less  than  2c.  In  such  a  situation,  if  the  first  failure  occurs 
after  a  long  enough  period  the  fee  pool  may  exceed  c  at  the  time  of  failure. 
Then  one  must  postulate  what  is  to  be  done:  the  contractor  may  take  up  to  c 
and  return  the  extra  amount  to  the  pool,  or  lose  the  extra  amount,  or  take 
the  extra  amount  as  a  sort  of  milestone  profit.  It  is  likely  that  a  study  of 
these  alternatives  would  produce  information  useful  in  the  structuring  of  actual 
incentive  warranties,  they  can  be  made  even  more  realistic  (and  involved)  by 
including  considerations  about  the  value  of  money.  However  for  the  present 
research  the  simplest  kind  of  model  was  used  in  which  the  contractor  receives 
the  fee  pool  amount  specified  by  the  incentive  fee  function  at  the  time  of  failure. 
At  that  time  the  fee  pool  starts  over  at  zero  and  increases  at  the  same  rate. 

For  some  contract  parameter  values  the  payment  could  go  negative  in  such  a 
model,  representing  a  profit  for  the  contractor.  This  situation  could  motivate 
the  contractor  to  assume  the  warranty  and  to  produce  equipment  with  a  desir¬ 
able  availability  profile.  However  such  an  instrument  may  also  improperly 
motivate  the  contractor,  particularly  if  parameter  values  are  negotiated  so  as 
to  give  an  advantage  to  the  contractor.  Thus  the  present  research  deals  with 
a  rather  simple  mathematical  model  form  that  contains,  by  implication  some 
important  contractual  considerations. 

This  incentive  fee  function  is  assumed  to  be  a  linear  function  having 
constant  slope  b.  If  Tj  and  T2  represent  the  random  variables:  time  to 
first  and  second  failures,  respectively,  then  the  payment  process  for  at  most 
two  failures  has  the  form: 

A2(t)  =  c-bT1  +  c  -  b(T2-  Tx). 
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Alternatively  if  R ^  and  R2  =  T 2  ■  T ^  are  the  random  failure  times 
to  the  first  failure  and  from  the  first  to  the  second  failure,  respectively,  the 


process  is  expressed  as: 


A2(t)  =  c-  bR1+c-b  R2- 


The  random  variables  R^  and  R2  are  commonly  used  in  availability 
work  and  will  form  the  two  random  variable  approach  used  in  much  of  this 


chapter. 

Let  f_  n  (r,,r,)  be  the  joint  probability  density  function  for  the 
R1R2  1  i 

random  variables  R^  and  R2  then: 

t  '-r2 

£[A2(t)]  =  J  f  [2c  -  +  r2)jfR^  (r1,r2)dr1  dr2 


t  '-r2 


E[A2(t)]  =  2c  p[Rl  +  R2<t]  -  b  Jf  (r1+  'jMr  R  (r  r  )(ir  dr 

OO 

A  similar  formulation  may  be  written  for  E[A2(t)]  from  which  the  Var[A2(t)] 
can  be  obtained  using  Var[A2(t)]  =  E[A2(t)]  -  (E[A;,(t)])^. 

In  the  rest  of  this  chapter  Rj^  and  R2  are  assumed  to  be  independent, 
exponentially  distributed  random  variables  so  that: 

,  ,  ,  1  >  .‘Vl  ”^2r2 

tR1R2,rl>r2)=Xl>'2e 


3.  2  Equal  Failure  Rates 


In  this  section  the  expected  value  and  variance  of  the  payment  process 
are  obtained  for  the  special  case  in  which  =  X2  =  X.  The  two  variable 
R  ,  R2,  representation  of  A2(t)  is  used.  For  convenience  in  calculation  let: 
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t  *-r2 


H(t)  =  f  f  (rx  +  r2)  fjtrj)  f2( r2)  drx  dr. 


o  o 


-^lrl  -X  r 

where  the  marginal  distributions  are  fj(rj)  =  e  ,  and  =  ^2e  ^ 


Thus: 


Direct  calculation  yields: 

p[ri  +  R2  <  t]  =  1  -  (1  +  Xt)  e"Xt 

H(t)  =  z[x  -  (x  +  1  +  T")  e‘Xt] 

E[A2(t)]  =  2cp[ri  +  R2<tJ  -  b  H(t) 


E[A2(';)]  =  2c  -  [2b(^  +  t  +^r>)  -2c(l+  Xt)]e’U  . 

If  the  expected  value  is  taken  for  the  full  range  of  t  (the  positive  real  axis) 

2b 

the  value  is  2  c  -  >y-  which  agrees  with  the  above  expression  in  the  limit  as 
t  —<o . 

2 

The  calculation  of  E[A2(t)]  is  more  involved  and  proceeds  as  follows: 

ECA^t)]  =  4  c2  p[rl  +  R2  <  t]  -  4bc  H(t)  +  b2  k(t) 
where  p("r^  +  *^2  —  ant*  **(t)  are  S^ven  above.  The  quantity 


t  t-r2 


k(t)  =  J  J  (rf  +  2rx  r2  +  r \)  X 


2.  2  ‘X(ri +  r2^ 


drl  dr2 


o  o 


Evaluation  of  this  integral  yields: 


k(t) 


■>LVCVf**2  ♦*£)•-“] 


-x6  'X 
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Therefore: 


EtA^t)]  =  4c2[l  -  (l+Xt)e'U]  -  8bc  -  Q-  +  t  +i^-)e‘Xtj 

.  ,,  2f  2  (  2  .  2t  .  .2  .  Xt3\  -Xtl 
+  3b  V^T  +  T  +  t  +^-Je  j 

Using  this  value  together  with  (ECA^t)])2  gives: 

(4c2X  -  8 b c  +l^)t 

+  b2t2  -  b2X  t3Je‘Xt  +  [.  2c(l  +  Xt) 

Xt2^  -2Xt 
+  2b^  +  t+-^-JI  e 

As  can  be  seen  this  rather  complicated  expression  has  a  constant  term  equal 
to  the  overall  variance  (obtained  here  as  the  limit  for  t -»<*>)  and  two  different 
exponential  terms.  The  later  terms  have  polynomials  in  t  as  coefficients. 
Detailed  study  of  the  expressions  for  E[A2(t)]  and  Var[A2(t)]  require 
numerical  approaches  outside  the  scope  of  the  present  study.  Their  com¬ 
plexity  indicates  that  a  corresponding  complexity  is  to  be  expected  in  the 
general  case  where  X^  i  X2*  The  Pre*®nt  special  case  serves  as  a  check  on 
the  more  general  expressions  developed  in  the  next  section. 

3,  3  Unequal  Failure  Rates 

As  indicated  by  the  equal  failure  rates  case  considered  above  the 
expressions  for  E[A2(t)]  and  particularly  for  Var[A2(t)]  become  complicated 
in  the  more  general  case.  For  this  reason  two  different  approaches  are  used, 
one,  the  single  variable  approach  being  simpler  to  carry  out  than  the  other 
which  is  the  standard  two  variable  treatment  of  the  kind  used  in  the  calculation 
of  the  previous  section. 
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Due  to  the  complexity  of  calculations  only  the  quantity  E[A.,(t)]  is 
obtained  by  the  two  variable  approach.  Its  value  serves  as  a  check  on  the 
calculations  using  the  one  variable  approach  and  also  for  consistency  with 
the  one  failure  result  E[A^(t)]. 

I  r»  Fy 

The  density  functions  for  and  are  X^e  and  \.,e  6  c 

respectively,  and 

E[A2(t)]  =  2c  P[R1  +  R2  <t]  -  bHz(t) 

t  t_r2  _x  r  x  r 

where  H2(t)  =  J  J  (rx  +  r£)  X^2  e  1  1  e  2  2  dru  dr2. 


t  t_r2 


+  S^*'Sri«'X2t2 


e  drj  dr2> 


Direct  calculation  yields: 


p[Rl  +  R2<t]  =  1  - 


-X2t  \2  .  -X2t  *Xjt 

6  +  xTrri  (e  -e 


n- 


X2t+^-e-xr+-^V(,+T7»fe ' 2  Xl>t-‘] 


These  values  give: 


:U2(t)]  =  2c  -  b(^+^)+J[b(t+^-+^)-2c]e  2 


2cX2  .  -Xgt  - X , t\  bX 


2  ( 

+*F^‘ 


•  (*♦£)( 


■  X-t  -  X .  t> 


e  -  e 


2  *i  r 


This  expression  can  be  written  in  several  ways  but  all  of  them  seem  to  be  of 
about  the  same  complexity.  The  direct  calculations  required,  though  simple, 
are  involved  and  it  does  not  seem  desirable  to  obtain  E[A^(t)]  and  hence  the 
variance  by  the  two  variable  method. 

Due  to  the  complexity  of  the  results  it  is  desirable  to  check  them 
against  each  other.  This  is  done  in  Section  3.4.  However  it  may  be  noted 
here  that  the  overall  expected  values,  obtained  if  the  process  goes  beyong  T 
should  agree  with  the  formulas  as  t-*«.  Let  =  2c  -b(R^  +  R2)  denote  the 
random  variable  when  time  is  not  limited  by  T  .  .Then 
E[A2]  =  2c  -  b^CRj]  +  E[R2]) 

E[A2]  =  2c  -  b()“  +  r^)  » 


and  Var[A2]  =  b^fvar[Rj3  +  Var[R2]^  , 


Var[A2]  =  b' 


These  values  agree  with  the  limit  values  for  the  Xj  =  X2  =  ^  case  and  may  be 
used  as  approximations  to  the  actual  expressions  for  E[A2(t)]  and  Var[A2(t)]. 

The  value  of  such  approximations  depends  in  part  on  the  probability  that  the 

4e 

failures  occur  after  T  .  Because  of  the  possibility  of  one  or  two  failures 
occurring  the  problem  of  estimating  the  degree  of  approximation  does  not  . 
appear  to  be  simple. 

3.  3.  2  The  One  Random  Variable  Approach 
An  alternative,  and  in  some  ways  simpler,  approach  to  the  study  of 
the  process  A2(t)  is  in  terms  of  the  single  random  variable  T2,  the  time 
(starting  from  aero)  to  the  second  failure  event.  The  reason  for  including 
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both  approaches  in  at  least  part  of  this  report  is  to  provide  optional  model 
forms  and  to  allow  calculations  based  on  one  approach  to  be  checked  against 
the  other  for  consistency. 

The  process  A2<t)  =  c  -  bT^  +  c  -  MT2  -  T^)  can  be  written  more 
simply  as: 

A2(t)  =  2c  -  hT2. 


To  study  A^(t)  in  terms  of  the  random  variable  T2  the  probability  density 

function  f^.  (v)  is  required.  Since  T2  =  R^  +  R2  the  previously  computed 
2 

value  for  P[R^  +  R2  <  v],  with  t  can  be  used  to  obtain  the  distribution 
function: 


Ft  (v) 
l2 


=  p[ri*R2<v]=1-.'1'2'' 


•X2V 


Differentiation  and  arrangement  of  terms  results  in: 


fT  (v) 
l2 


XlX2  {  ~  Xlv  -x2*\ 


Partial  expectations  of  T2  restricting  v  <  t  are  used  to  obtain  expectations 
for  A~(t).  The  concept  of  partial  expectations  is  discussed  more  fully  in 
Chapter  4. 

Let  Et[T2]  denote  the  partial  expectation  of  T2,  restricting  v  to: 
v<t.  Then: 


W&j<' 


V 


-  x.x. 
E«[T2]  =T-T7 


V 


17 
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Then  E[A2(t)]  =  2cP[T2<t]  -bE^[T2]/  which  gives: 


E[A,(t)]  =  2c 


1  -e 


-X2t  X 


+ctt( 


-X2l  -Ll 
e  -  e 


2  "1 


bXlX2|  1  1  J 
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-X-t  -X,t  -X,t  -X.t 
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■1 


2  *1 

This  expression  can  be  shown  to  equal  the  expression  obtained  in  Section  3.  3.  1 
and  in  the  limit  X2-*X^  =X  to  yield  the  expression  from  Section  3.  2,  these 
considerations  are  discussed  in  Section  3.4. 


Z  2 

To  obtain  E[A2(t)]  and  hence  Var[A2(t)]  the  expression  E  [T^]  is 


required. 


Trr^2n_XlX2  r  2'  "Xlv  ”X2VX , 

Et[T2]  -  rpx7' v  -e  )dv 


v  rr2l  X1X2 

e.Ct2]  -1757 
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ECA2(t)]  =  Et[4c2  -  4bcT2  +  b2  T2] 


=  4  c2  PCT2  <  t]  -  4bcEt[T2]  +  b2  Et[T2] 


All  the  quantities  required  for  this  expression  have  been  obtained  above.  By 
putting  them  into  the  expression  it  becomes: 
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where 


E[A£(t)]  =  D{  +  D2  +  D3  , 


Dx  =  4c 
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These  expressions  can  be  used  to  give  the  variance  in  the  form: 

Var[A2(t)]  =  E[A2(t)]  -  (eCA.,^)])2  . 

The  only  way  to  consider  such  complex  expressions  is  as  numerical  functions 
of  t.  This  kind  of  detailed  calculations  is  outside  the  scope  of  this  report. 
However  internal  checks  and  limit  cases  can  reasonably  be  considered  as  is 
done  in  the  next  section.  However  the  case  where  t-*«  can  be  caluclated 
directly.  This  yields: 


E[A2]  =  2c  -  b(^-  + 


2 

—  .2,  .  2  X1X2  (  1  1  .  b  X1X2  /  2  Z\ 

EtA^-4c 


=  4  c  -  4  be 


(^+^)+2”2^+T7i+i|-) 
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These  values  give: 


2f  1  i  "> 

Var[A2]  =  b  | — j  which  agrees  with  the  overall  expected  value 

~X ,  X  - 


1  'v2 

and  variance  values  obtained  at  the  end  of  Section  3.  3.  1. 


3.4  Verifications  and  Limit  Values 
3.4.  1  Agreement  of  Approaches 

In  Section  3.3.  1  the  expression  for  E[A2(t)],  denoted  as  E1[A2<t)] 


was: 


The  expression  obtained  in  Section  3.  3.  2,  denoted  as  E2[A2(t)]  was: 


E2[A2(t)]  =  2c]  1  -  e 


-X2t 


-X,t 


-X.t 


X2"X1 


-  e 
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bV2 


X2'^l 


1 

,  /  -v  > 

i  ,  t(®  e 

-x2t 

"X  ^ 

x  e 
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>7 

x  x2  \l 

rv»  oj 
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X  ^ 
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It  is  not  immediately  evident  that  these  two  expressions  are  the  same.  Showing 
that  they  are  the  same  gives  an  indication  that  the  expression  for  E[A2<t)]  is 
correct  and  that  the  two  approaches  to  model  formulation  and  analysis  are 
equivalent.  This  constitutes  an  internal  consistency  check  on  the  research. 

To  show  that  the  two  expressions  are  the  same,  similar  terms  are 
compared.  The  constant  terms  are  both  equal  to: 

2c  -  b^-+^-^  .  The  remaining  term  with  c  as  a  factor  is  also  seen 

to  be  the  same  in  both  cases,  equal  to: 
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The  remaining  term  in  E^CA^t)]  is: 


In  E^CA^t)]  this  term  corresponds  to: 


~b  ^1^2 


A  little  algebra  shows  that  both  of  these  expressions  are  equal  to: 


Therefore  it  is  established  that  E^[A2(t)]  =  E2[A2(t)]  as  required  for  con¬ 
sistency.  The  two  approaches  give  the  same  value  for  E[A2(t)].  Since  the 
variance  was  only  obtained  for  the  one  variable  approach  no  check  cf  tnis 
kind  can  be  made  for  its  expression. 

3.4.2  Limits  as  Failure  Rates  Become  Equal 
Both  E[A-,(t)]  and  Var[A2(t)]  are  available  in  the  case  where 
=  X2  =  X,  from  Section  3.2.  The  expressions  for  these  quantities  for  the 
general  case  X ^  t  X2  are  obtained  in  Section  3.  3.  2.  These  should  tend  to  the 
special  case  values  as  X^  -  X2  tends  to  zero.  The  expressions  in  question 
are  simple  functions  of  the  failure  rates  X^  and  X2  and  the  limit  as  Xj  -  X2  -»0 
certainly  exists.  The  intent  of  this  section  is  to  evaluate  that  limit  and  show 
that  it  gives  the  desired  special  case  value.  To  do  this  it  is  satisfactory  to 
take  any  particular  approach  of  the  variables  to  their  limit  values.  Two  forms 
are  used  in  the  following:  X2  -  X^  =  v-’®  and  ^2  wb^e 
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I 


First  consider  the  E[A£(t)]  expression.  Write  this  as; 

=  2cP  -  bF^t).  The  limits  of  each  part  will  be  obtained  and  shown 
to  equal  the  corresponding  P  and  H(t)  expressions  in  the  Xj  =  case. 


P  = 


-X  t  . 

e  4-  X 


-Xt 


lim 

VO 


-  1 


P  =  1  -  (1  -  Xt)  e”Xt  t  which  is  the  special  case  expression. 


Note  that 


lim 

Vh-° 


1  +  e-(x2-xi>* 
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=  -  t  so  that 


lim 


(‘ +  r)  e'Xl  +  x  e"u(t  +  *•)<-*> 


2 

X 


+  t  + 


the  special  case  value. 


Therefore 


lim 

VX2~° 


EtA^Jt)]  =  the  special  case  value  of  EtA^t)]  when  X^  =  X2=X  . 


Similar  analysis  is  used  to  show  that  the  general  case  variance  goes  to 
the  special  form.  The  procedure  is  to  show  that  the  expression  for  E[A2{t)j 
has  the  correct  limit  form  since  it  has  already  been  established  that  the 
expected  value  has  the  proper  limit.  The  rather  involved  calculations  are 
carried  out  for  each  of  the  terms  D^,  D£»  and  used  to  express  E[A2(t)] 
ia  Section  3.3.2.  These  may  be  compared  individually  with  the  corresponding 
terms  in  the  special  case  expression  for  EfA^t)]  given  in  Section  3.2.  Direct 
calculation  of  limits,  setting  X2  =  X  and  letting  Xj-»X  yields; 

lim  D,  =  4c2  j  1  -  e"Xt  -  Xte*Xtj,  which  ia  the  correct  equal  failure 

Xj-*X  1  L 


-TT 
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rate  term  containing  c  . 

lim  =  -  4b  c!^  -  (~  +  2t  +  X  t^  j  e  ^ I ,  which  is  the  correct  term 
X  ^  -»X  ^  '  -* 

containing  be. 

p.  ,  2  r 6  / ,.2  ,  ,  .3  .  6t  .  6  A 

Xj-X  3  l X L  '  x  x£/ 

2 

the  correct  b  term.  Thus  the  VarCA^t)]  for  the  general  failure  rate  case 
does  become  the  value  obtained  for  the  equal  case  X^  ~  ^2  =  X.  This  provides 
a  further  consistency  check  on  the  models  and  their  analysis. 

3.4.  3  Limits  as  the  Second  Failure  Rate  Increases 
If  X^  becomes  large  the  time  between  the  first  and  second  failures 
reduces  in  a  probabilistic  sense.  In  the  limit  as  X^-*®  the  two  failures  occur 
together.  In  this  case  the  results  should  agree  with  the  one  failure  model 
except  that  the  cost  parameter  becomes  2c  corresponding  to  the  fact  that 
there  are  actually  two  failures,  or  double  the  cost,  at  the  combined  failure 
time.  Evaluation  of  this  limit  case  serves  as  a  check  on  the  single  and  double 
failure  models. 

Section  2.3  gives  the  following  values  for  the  one  failure  case.  Para¬ 
meters  are  adjusted  so  as  to  correspond  to  the  model  of  the  present  Chapter 
(aQ  =  0,  t'  =  0,  and  c  =  2c). 

ECAjIt)  .2c(l-  .-Xt)  -  b[i  -  (£  +  t)e-Xt] 


V.^AjW]  -£)  -b2 


.2  b2  -\t 

*  T*J* 


.(zc-i-bt)  . 


-2Xt 
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2 

Using  the  expressions  for  E[A2(t)]  and  E[A2(t)]  developed  in  Section  3.3.2, 
taking  the  limit  as  X2“*®  yields  the  following 

lim  E[A2(t)]  =  2c[l  -  e  XtJ  -  bX ~\~%  +  yj  e  ^  | 

x2—  Lx  x  J 

where  Xj  =  X  for  comparison  with  the  one  failure  expressions.  So  that 

lim  E[A2(t)]  =  ECA^t)]. 

X2  — 

In  the  two  failure  case  set  X^  =  X  and  let  X2—®  to  obtain: 

..  rrA2mi  a  2!~i  -Xt]  ..  *  1  t  -Xt  e~Xt 
lim  E[A2(t)]  =  4c  |_1  -  e  J  -  4 be X  — j  -  »j— 

X2"®  X  X 


+  b2x  -  y-  e'Xt  +  -  2(t  +  i)  2-z- 

m 


For  the  one  failure  case  the  expression  (with  c  replaced  by  2c)  is: 

E[Aj(t)]  =4c2  -i^+-^.-^4c2-4bc(tt^)  +  b2(t*^) 

2  2 

In  this  form  direct  comparison  of  the  c  ,  be,  and  b  terms  show  that: 
lim  E[A2<t)]  =  E[A2(t)]  . 

x2  — 

It  follows  from  these  results  that  lim  Var[A2(t)]  =  Var[A^(t)]  as  desired. 

x2-» 

Thus  the  limit  connection  between  the  one  failure  and  two  failure  cases  is 
established. 
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3.  5  Disc  us  s  ion 


The  expressions  developed  in  this  chapter  represent  a  basic  mathemat¬ 
ical  methodology  for  the  study  of  incentive  warranties  when  there  are  at  most 
two  failures.  The  results  are  complicated  functions  of  the  contract  parameters, 
failure  rates,  and  time.  They  can  best  be  studied  in  detail  by  numerical  evalua¬ 
tion  and  graphing  of  the  representational  quantities  E[A2<t)],  and  Var[A2<t)]. 
However  special  studies  can  be  made  without  detailed  numerical  calculation. 

Such  studies  can  address  some  specific  questions  of  interest  and  provide  in¬ 
sight  about  the  model.  A  typical  special  consideration  is  given  in  this  section, 
both  for  its  intrinsic  interest  and  as  an  example  of  the  kind  of  study  that  can  be 
made.  Some  considerations  about  the  present  value  of  money  in  the  more  than 
one  failure  case  are  also  discussed.  Some  additional  checks  on  complicated 
expressions  using  dimensional  analysis  are  discussed  at  the  end  of  this  section. 


In  Section  3.  3.  2  an  expression  is  given  for  Et[T2],  the  expected  value 
of  the  time  to  second  failure  up  to  time  t  when  Xj  t  \^.  If  this  expression  is 
considered  for  the  limit  case  with  \  and  X^  —X  one  obtains: 


E.[T,]  =  £  +  *  11x11  — 


-Xt  ~^1* 
X,  e  U-Xe  1 


Xj-X 


X  -  X, 


+  lim 


,  2  -Xt  ,2 
Xj  e  -  X  e 


X.  -X  X,  (X  -  X.) 


EtCT2]=f-[2t  +  u2+f]e' 


A  special  case  of  some  interest  is  obtained  when  the  expected  time  to 
failure  is  T  / 2  so  that  in  a  purely  expected  value  sense  the  first  failure  occurs 
half  way  through  the  warranty  period  and  the  second  failure  occurs  at  the  end. 
This  can  be  thought  of  somewhat  like  the  "one  hoss  shay"  case  of  Chapter  2, 
with  of  course  some  major  differences  due  to  the  exponential  nature  of  the 
failure  process.  For  this  case  X  =  2/T*  and  the  expected  payment  by  the  end 
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of  the  warranty  period  is  obtained  by  setting  t  =  T  in  the  appropriate  expres¬ 
sions.  For  this  purpose  one  could  use  the  quantity  E[A2(t)]  with  X  ^  =  X2  ~ X* 
However  the  same  result  can  be  obtained  by  the  following  alternative  pro¬ 
cedure: 

E[A2(T*)]  =  2c  P[T2  <  T*]  -  bE^[T2]  , 

Jjc 

since  P[T2<  T*]  =  1-3  e'XT  ,  and  E  *[T2]  =  T*  -  5T*e'XT  ,  the 

result  becomes  E[A2(T*)]  =  2c  -  b  T*  +  [5  T*  -  6c]  e 

When  the  second  failure  occurs  at  T  the  payment  is  equal  to 

jJ; 

2c  -  b  T  no  matter  when  the  first  failure  occurs.  Thus  if  there  is  a  deter- 

$ 

minjstic  second  failure  the  total  payment  will  be  exactly  2c  -  b  T  .  The 

* 

expected  payment  when  Xj=X2  =  2/T  will  be  greater  or  less  than  the 

$ 

deterministic  amount,  depending  on  whether  5  b  T  -  6c  is  positive  or  negative 
respectively.  This  observation  provides  an  opportunity  to  relate  the  contract 
parameters  b,  c,  and  T  .  However  the  relation  only  holds  for  the  special 
failure  profile  governed  by  the  common  failure  rate  2/T  . 

By  having  the  second  failure  occur  at  T  ,  given  that  two  failures  must 
occur  in  the  warranty  period,  the  contractor  achieves  maximum  benefit  from 
the  fee  pool.  Therefore  he  will  expect  to  pay  no  less  than  the  amount 
2c  -  b  T  .  Hence  one  obtains  5b  T  -  6c  >  0.  It  is  a  basic  assumption  of  the 
present  model  methodology  that  the  contractor  can  not  profit  from  the  failure 
events,  imposing  the  condition  2c  -  b  T*  >  0.  These  two  conditions  lead  to  a 
rather  tight  relation  between  the  three  contract  parameters,  of  the  form: 

2c  >  b  T*  >  1.  2c  . 

* 

If  the  contractor  is  willing  to  have  the  expected  payment  less  than  2c  -  b  T 


he  will  have  to  assume  a  greater  risk  of  having  to  pay  more  than  the  expected 
value.  This  is  a  less  conservative  position  than  the  one  described  above.  It 
imposes  the  condition:  5b  T  -  6c  <  0  which  reduces  to:  b  T  <  1.  2c  which  is 
always  consistant  with  the  nonprofit  assumption  b  T  <  2c. 

This  kind  of  analysis  can  serve  as  a  guide  for  the  selection  of  contract 
parameters,  with  proper  caution  for  the  underlying  assumptions  about  the 
random  failure  profile.  The  cost  c  and  warranty  period  T  are  relatively 
easy  to  establish  so  that  a  condition  such  as  2c  >  b  T  >  1.  2c  provides  a  means 
of  determining  the  incentive  pool  slope  b  by  satisfying: 

T  T 


When  the  present  value  of  money  is  considered  a  number  of  questions 
arise.  For  example  if  the  first  failure  event  occurs  at  a  time  for  which 
c  -  b  Tj  <  0  there  is  an  excess  pool  value  which  might  be  invested  in  such  a 
way  as  to  provide  more  future  value  than  simply  retaining  it  as  part  of  the 
pool  (the  procedure  used  in  the  models  of  this  chapter).  Such  excess  pool 
investment  would  have  to  be  allowed  or  disallowed  by  the  warranty  contract. 
One  consideration  in  such  a  situation  would  be  the  possibility  of  allowing  the 
excess  as  an  award  fee  in  the  event  of  no  future  failure,  or  of  requiring  it 
to  be  used  only  when  and  if  a  second  failure  occurs.  These  are  interesting 
possibilities  for  motivational  contracting  in  these  situations.  Another  possi¬ 
bility  is  for  the  contract  to  allow  the  contractor  an  option  as  to  how  much  of 
the  incentive  fee  pool  he  uses  toward  payment  of  the  first  failure.  In  a  dis¬ 
counting  (inflationary)  environment  one  can  argue  that  the  full  amount  available 
should  be  spent  at  the  failure  time  since  it  will  only  be  worth  less  at  any  future 
expenditure.  However  if  the  contractor  is  allowed  to  invest  a  part  of  the  fee 
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pool  that  is  due  to  him  at  the  failure  time  then  he  can  balance  investment 
income  against  the  discount  rate  effect,  and  by  assuming  various  levels  of 
risk,  hope  to  improve  his  position  for  payment  of  the  second  failure  if  it 
occurs.  When  no  second  failure  occurs  the  contract  must  specify  the  dispo¬ 
sition  of  such  funds.  Still  another  possibility  is  to  allow  the  contractor  to 
borrow  against  future  levels  of  the  fee  pool.  This  is  one  way  to  overcome 
discounting  effects  by  spending  money  at  an  earlier  time.  There  are  problems 
with  this  concept  such  as  how  to  assess  the  contractor  in  the  case  of  no  second 
failure.  A  possibility  seems  to  arise  which  may  provide  a  negative  motivation 
in  which  the  contractor  will  actually  desire  a  second  failure  in  order  to  release 
the  fee  pool  funds.  More  detailed  investigation  of  value  of  money  concepts  is 
beyond  the  scope  of  this  report.  However  they  may  be  expected  to  play  an 
important  role  in  actual  warranty  contracts  and  lead  to  complications  requiring 
careful  analysis  by  both  contractor  and  customer. 

The  expressions  developed  in  this  chapter  are  complicated  and  it  is 
important  to  check  them  as  thoroughly  as  possible  against  computational  or 
theoretical  error.  A  major  class  of  such  checks  is  given  in  Section  3.4.  An 
additional  type  of  check  is  to  consider  the  dimension  of  each  term  in  an  ex¬ 
pression  to  be  sure  they  all  represent  the  same  quantity.  This  kind  of  dimen¬ 
sional  analysis  check  is  illustrated  below  for  E[A2(t)]  in  the  equal  failure 
rate  case.  The  quantity  to  be  checked  is: 

E[A2(t)]  =  2c  +^2b(i  +  t+^-)  -  2c (1  +Xt)je'Xt 

The  dimensions  of  the  various  quantities  are  given  in  the  following  table 
where  T  denotes  "time"  and  Y  denotes  dollars: 
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Quantity 

c 

b 

\ 

t 

Dimension 

Y 

Y/T 

1/T 

T 

Since  \t  is  dimensionless  the  factor  e  ^  is  dimensionally  correct.  The 
expression  is  now  put  into  "dimensional"  terms: 


dimension  E[A,(t)] 


=  Y 


» 

11+  I 

T  T 


+  T 


+  ! 
T 


_2\ 
T  l 


-  y(  1  + 


It)  . 

T  / 


It  is  seen  that  every  term  has  the  dimension  of  dollars  as  should  be  the  case. 
This  kind  of  check  is  particularly  useful  in  finding  mistakes  such  as  missing 
factors  or  incorrect  powers  of  variables. 


3.20 


Chapter  4.  General  Availability  Case 

Section  4.  1  General  Methodological  Model 


4.  1.  1  Model  Form 

In  this  chapter  the  payment  process  A^(t)  for  k  failures  is  considered 
together  with  the  general  payment  process  A(t)  which  does  not  depend  on  any 
specific  assumption  about  the  number  of  failures.  It  is  important  to  observe 
that  the  payment  processes  are  considered  up  to  a  time  t  which  is  of  most 

3$C  3$C 

interest  when  t  <  T  ,  with  T  as  the  warranty  period.  The  expression  of 
results  as  a  function  of  time  t  gives  considerable  structure  to  the  availability 
warranty  models  and  is  imposed,  in  any  case,  by  the  finite  warranty  period 
T  .  However  such  a  model  representation  causes  restrictions  in  both  the 
results  and  analysis  techniques  normally  available  for  the  study  of  failure 
processes.  In  the  usual  case  there  is  no  restriction  on  the  time  and  full  dis¬ 
tributions  of  failure  times  can  be  employed.  This  is  not  the  case  in  the  present 
study  as  already  observed  in  previous  chapters  for  the  one  and  two  failure 
situations.  The  effect  of  the  more  complicated  model  structure  imposed  by  the 
finite  time  condition  is  even  more  evident  in  the  general  considerations  of  this 
chapter. 

The  general  methodological  model  for  incentive  availability  warranties 


utilizes  the  following  quantities: 

F^.  (t)  =  P^T^  <  t j,  the  distribution  function  of  the  k^,  failure,  denoted 

by  the  random  variable  T^. 

f_  (v)  =  the  density  function  for  T. 

Tk  * 

eJ"TjJ  =  the  partial  expectation  of  T^  restricted  to  the  range  T^  <  t: 
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t 

Et':TkJ=i  vfT Jv)dv 
0  k 

P 

E!  A^(t)j  =  expected  payment  under  the  warranty,  by  time  t  <  T  , 
when  there  are  at  most  k  failures. 

Ak(t)  the  payment  process  up  to  time  t  when  there  are  at  most  k 
failures. 

One  could  of  course  also  define  the  variance  of  the  payment  process, 
Var|~Ak(t)J,  but  for  k>2  this  quantity  is  an  extremely  complicated  function 
of  t  and  the  contract  parameters.  Therefore  in  the  general  case  it  is  not 
employed  as  a  representational  feature  of  the  model. 

A(t)  is  the  payment  process  up  to  time  t  no  matter  how  many  failures 
may  occur. 

E[A(t)]  is  the  expected  value  of  A(t)  up  to  time  t 

N  =  the  number  of  failures  that  occur 

P[N  =  k]  is  the  probability  function  for  the  discrete  random  variable  N. 

In  the  general  methodology  the  random  variable  N  is  defined  indepen¬ 
dently  from  the  failure  time  random  variables  T^.  Of  course  the  failure  pro¬ 
cess  imposes  a  relation  between  the  failure  times  and  the  number  of  failures. 
For  example  the  well  known  generation  of  exponential  failure  times  by  a 
Poisson  failure  process.  However  in  the  present  methodology  the  point  of  view 
is  different.  As  a  contract  negotiation  decision  position  the  parties  to  the  con¬ 
tract  may  assume  some  a  priori  probability  function  P[N  =k].  The  generation 
of  failures  by  the  failure  time  process  is  unrelated  to  such  an  a  priori  speci¬ 
fication  of  number  of  failures  for  the  warranty  period  (or  indeed  for  any  time 
value  t).  In  fact  the  generation  process  and  the  specification  of  failures  may 
be  related  in  various  ways  that  involve  considerations  outside  the  scope  of  the 
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esearch  reported  on  in  this  paper  but  falling  within  the  considerations  of 
future  research  as  mentioned  in  Chapter  5. 

The  above  definitions  and  considerations  lead  directly  to  the  formula¬ 
tion  of  the  general  model  in  the  following  form: 


E[A(t)]  =  l  E^t)]  P[N  =  k]  , 
k=l 

E[Afc(t)]  =  kcpfTk<t]  -  b  Et[Tk]  , 


t 

r 

o 


f_,  (v)dv  . 
ik 


In  these  models  E[A{t)]  and  eTAj^JJ  are  partial  expectations  taken 
over  the  range  [0,t],  generalizing  similar  results  from  Chapters  2  and  3. 

If  R.  is  the  random  time  from  the  i-1  st.  to  the  ith  failure: 

Tk  =  +  R2  +  .  . .  +  Rk  and  f^,  (v)  can  be  generated  from  the  individual 

failure  time  distributions  f_  (v).  In  this  report  the  only  specific  failure  time 

1 

random  variables  considered  are  exponential. 

4.  1.  2  Approaches  to  Analysis 

Because  Tk  is  the  sum  of  k  independent  random  variables  one  might 
consider  the  possibility  of  computing  the  moments  for  Tk  (or  indeed  its 
distribution  form)  by  using  characteristic  (or  moment)  generating  functions. 
Equivalently  sums  of  expectation  formulas  suggest  themselves.  However  these 
classical  approaches,  so  useful  in  the  usual  full  expectation  situation  become 
so  involved  in  the  partial  expectation  situation  as  to  be  of  little  value.  There¬ 
fore  it  is  necessary  to  proceed  directly  to  the  calculation  of  desired  quantities. 
Before  proceeding  to  consideration  of  such  calculations  some  results  will  be 


‘.7-  '>•««,  «.n»3V<  -s 
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given  for  partial  expectations  to  illustrate  the  problems  they  present  for 
analysis  and  to  record  their  form  as  a  by-product  of  the  present  research. 

It  should  be  explicitly  recognized  that  there  is  nothing  technically  wrong  with 
using  the  proper  form  of  partial  expectation  quantities.  The  problem  is  not 
theoretical  but  practical  application  of  the  complicated  forms  involved. 

As  the  first  example  of  partial  expectation  calculations  consider  the 
random  variable  T^,  time  to  second  failure.  Since  T£  =  Rj  +  R2  full 
expectation  follows  the  well  known  formula  E^T2j  =  EpR^j  +  £^2^  in  which 
the  marginal  distributions  for  R^  and  R2  can  be  used  to  compute  E^R^J  and 
E^R2  .  In  terms  of  partial  expectations  a  formal  similarity  exists  as  may 
be  seen  from  the  following  development: 

t  t-r2 

Et[T2W  I  'rl  +  r2>£R..B  <rl-r2,drI  dr2 
00  1  L 


t  t_r2 


t  fc-r2 


=  JJ  rlfR.,R,(rl,r2)drl  dr2  +  ,f  I  r2  fR.,R2(rl'r2)drl  dr2 
0  0  1  L  0  0 

*  *t[*i] + "tM  • 


The  derivation  of  this  formal  result  clearly  indicates  the  need  to  compute 
Ef^RiJ  in  the  mvilti variate  sense  due  to  the  partial  expectation  effect.  This 
considerably  complicates  the  use  of  these  formulas. 

Though  the  use  of  partial  expectations  complicates  calculations  they  do 
give  correct  results  when  properly  employed.  As  an  illustration  in  the  expo¬ 
nential  failure  case  with  Rx  and  R2  independent  random  variables  having 
failure  rate  parameters  and  \2  respectively,  direct  calculation  yields: 
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In  Section  3.3.2  the  following  result  is  given: 


A  little  algebra  shows  that  these  two  expressions  are  the  same.  The  example 
illustrates  the  fact  that  the  classical  sum  of  expectation  procedures  are  not 
particularly  helpful  in  dealing  with  partial  expectations . 

Because  of  the  increasing  complexity  of  the  calculations  for  more  than 
two  failures  it  would  be  desirable  to  use  characteristic  function  methods  to 
compute  variance.  However  because  of  the  partial  expectations  these  methods, 
so  powerful  in  classical  application,  do  not  help  in  the  present  application. 

This  may  be  observed  from  the  following  three  considerations. 

k 

1.  The  usual  product  formula  cp  (z)  =  H  (z)  for  independent 

1 k  j=l  i 

random  variables  R.  does  not  hold  at  ail  for  the  partial  characteristic  function 

J 
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for  any  random  variable  X. 


defined  as  C0„(z,  t)  =  E.i  e*r^  ! 

A  t  c.  -j 

2.  The  composite  partial  characteristic  functions  :p_  (z,t)  are 

j 

complicated  to  compute  since  they  must  be  evaluated  in  the  full  multivariate 
environment,  i.  e. 


-X  r 


•(X.  -  iz)r. 


cp_  (z,t)  =  |  ...  i  II  \  e  n  ndr  X.  e  J  Jdr. 

Rs  :.*L °  _  ; J 

ntj 


^  r  2,  *  •  •  '  ^  ^  * 


e.g.  when  k=3 


t  *-'1  t_rl"r2 


i»  r 

(*»*)  a  J  j 

1  0  0  0 


J' 


X2r2  _X3r3  ,  .  '(X  i-izjr, 

X2X2  e  e  dr^dr,  X,  e  dr,, 


2  3  1 


1’ 


which  even  for  this  relatively  simple  case  (k=3)  is  not  at  all  simple  to  evaluate. 

3.  The  useful  formulas  for  expectation  and  most  particularly  for 
variance  no  longer  hold  in  the  partial  expectation  situation. 

t 

cpT(z,t)  =  Et|_elzT]  =  j  e1ZVfT(v)dv 

0 

2 

=  P[T<  t]  +  Et[T]  iz  +  Et[T2]  +  •  •  • 


so  that  cp^.{0,t)  =  P[T<t]  unlike  the  classical  characteristic  function  for  which 
cPx(O)  =  1  greatly  simplifying  calculations. 

The  characteristic  function  may  be  used  to  compute  the  variance  directly 
by  the  classical  result: 


d  log  =Px(z) 

1 - 


dz 


z=0 


=  -  Var  [X]  . 


For  partial  expectations  this  result  is  greatly  complicated.  Direct  calculation 


shows: 


Vart[T]  ={^j(  EtCT]) 


2  log=PT(z,t) 


-  P. 


dz 


z  =  0 


d  log  <P  (Z,t) 
where  iEt[T]  =  P - ^ - 


z=0 


and  P  =  P[T  <  t]  . 


The  above  considerations  of  partial  expectations  show  that  the  most 
reasonable  approach  to  analysis  of  the  general  failure  case  is  by  means  of 
direct  calculation  on  the  random  variables  T,_.  Because  of  the  complexity 
involved  evaluation  of  Var^T^J  and  hence  Var["A^(t)J  must  be  omitted. 

Tk 

k=3,  the  results  would  be  extremely  cumbersome.  This  is  true  even  for  the 
equal  failure  rates  case  which,  though  simpler  than  the  general  case,  is  still 
complicated  as  shown  in  Section  4.  3. 

The  analysis  will  form  the  density  function  f^.  (v)  for  the  random 
variables  T^  and  use  it  to  produce  the  required  quantities  P{_Tk<tJ,  E^TjJ  ’ 
and  E^A^t)!.  These  may  then  be  combined  with  the  a  priori  values  P[N=k] 
to  form  E[A(t)]  the  unconditional  warranty  payment  process. 

Section  4.  2  Exponential  Model  With  Unequal  Failure  Rates 
In  this  section  the  times  between  failures  R.  will  be  independent 
exponentially  distributed  random  variables  with  failure  rates  Xj-  The  time  to 
the  k**1  failure  is  given  by  T^  where  T^  =  Rj  +  &2  +  •  •  •  +  R^*  To  start  the 
analysis  the  random  variable  T^  is  considered.  Though  the  calculations 
involved  are  elementary  in  form  they  are  rather  complicated  in  execution  and 


Though  it  is  feasible  to  compute  Var^ 


for  some  specific  cases,  such  as 


only  the  major  steps  will  be  given. 

t  **1 


'  l"r2 


•T  (t) 
A  3 


’[T3i,J  *  '  J  J  Ws  * 


-X3r3-X2r2-X1ri 


dr3dr2  dr  x 


0  0 
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Direct  calculation  and  some  algebraic  arrangement  yields: 


~X , t  X  ,  X  ~ 

FT3(t)  =  1-e  + 


f  ”X2*  “Xi^  J  ^  1  X2  I  ■X'jt] 

\e  "e  r  ('X243)(x1A3Ae  *e  / 


This  distribution  function  for  T3  has  the  proper  limit  of  one  as  t  becomes 

unbounded.  The  density  function  for  is  obtained  as  the  derivative  of 

Ft  (t).  It  can  be  expressed  in  the  form: 

A3 


fT3(t)  =  xix2x3  (x1-x2y(x1-x3i +  (xp^xpqy +  Tpr 


As  a  check  on  this  result  it  may  be  observed  that: 


-x3t  1 

a 

iropqi 


^1^2  (  ”X2£  “ X  |  ("t)  , 

iim  fT  (t)  =  Ve  *  -  e  1  /=fT(t)  . 

X 3  “*oo  l3  xl~K2  l2 


Next  the  general  form  for  f_  (t)  will  be  developed.  It  may  be  observed 

Ak 


t  t-w 


^k^1] =  p[Tk-i +  Rk^] =  r  j  fTl  <w>  xk  *  k  dr  dw 


At  this  point  it  is  helpful  to  introduce  the  simplifying  notation  f_  (t)  =  f,  (t). 

ik  K 

Then  upon  partial  integration  the  equation  becomes: 


(•  r  “Xjc(t-w)-j 

FT  ^  =  J  Vl^L1  "  ®  -idw* 


Applying  Leibniz's  rule  to  differentiate  the  integral  expression  yields: 


f  -X.  (t-w) 

W)=\  £k-i(w)xke  dw 


which  determines  the  densities  f^(t)  as  solutions  to  an  integral  difference 


>:*■* 


equation.  The  solution  is  obtained  by  the  guess  and  verify  technique.  Thus 
it  may  be  verified  by  insertion  into  the  above  equation  that  the  solution  is  given 
by  the  following  expression: 


\  e 

yt)  =  xr..xk  i  (x1-\i)a2-xi)...(xk-xi) 


in  which  the  denominator  products  contain  k-1  factors,  with  the  zero  factor 
XrXi  absent.  This  result  generalizes  the  quantity  yt)  previously  obtained 
directly  as  the  density  function  for  T^. 

The  quantities  of  interest  for  the  incentive  warranty  model  may  now  be 
calculated  using  the  density  function  fk(t).  Using  the  definitions  of  the  quantities 
involved  and  simple  integration  yields: 


$  ( 


1  -  e 


-XA 


P[Tk-t]  =  Xl***Xk  L  X.(X1-Xi)...(Xk-X. 


k  ,  "Xi*  .  “Xifc 
Pi  ^  1  -  e  -  tXj  e 

EtLTkJ  =  xr  *-xk  L  TTT.  “  “  TT 


i=l  Xi  ^Xl"Xi^*  ’  *  ^Xk”Xi^ 


Since  there  are  k  failures  the  repair /replacement  cost  is  kc  so  that: 


E[Ak(.)]  =  kcp[Tk<t]-bEt[Tk] 


This  result  yields  the  previous  results  for  eTa ^(t)J  and  E^AjtttJ  from 
Chapters  2  and  3  respectively.  One  value  of  the  previous  alternative  derivations 
is  that  they  provide  checks  on  the  general  result  presented  here. 

See  complement  to  Chapter  4. 
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Section  4.  3  Exponential  Model  With  Equal  Failure  Rates 
In  the  context  of  warranty  modeling,  as  in  many  other  aspects  of 
availability,  the  special  case  of  equal  failure  rates  for  each  failure  is  likely 
to  be  the  most  realistic.  It  is,  of  course,  also  simpler  in  form  than  the 
general  failure  rate  model  of  Section  4.  2.  Intermediate  situations  in  which 
some,  but  not  all,  failure  rates  are  different  are  not  considered  in  this  report. 

When  =  Rj  +  ...  +  R^  and  R.  are  independent,  exponential  random 
variables  with  the  same  failure  rate  X  it  is  well  known  that  the  density  function 
for  is  a  gamma  density.  A  notation  specifying  the  equal  failure  rate  case 
is  obtained  by  placing  a  diacritical  mark  over  the  function  symbol. 

Even  though  the  density  function  is  known  to  be  gamma  it  may  be 
obtained  from  the  integral  difference  equation  in  the  same  way  the  general 
density  was  obtained  in  Section  4.  2.  From  results  in  Chapters  2  and  3  the 
values  for  and  T£  are  already  available  as: 

fjW  =  X  e_Xv 

?2M  =  X2  ve"^v  . 

The  integral  difference  equation  yields: 

f>)  =  x3  £  .-Xv  . 


The  result  for  arbitrary  k,  which  may  be  shown  to  satisfy  the  integral  dif 
ference  equation  is: 


X1^'1 

t 


e 


Xv 


In  the  equal  failure  rates  model  a  major  role  is  played  by  the  incomplete 
gamma  function  Y(a,  x)  defined  as: 
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X 


Y(a,x)  =  J  e  v  v3  *  dv  . 


The  quantities  required  for  the  model  may  now  be  developed. 
p  n  k  k- 1  *\v  j  i 

p[Tk -  -  J  fk.-nf — -  -  nrnrr  v<k-kt>  • 


p«tTkJ  *  J  ferr  •'Xw*»  *  t  ^[Tk+1 <  *] 

0 

®tLTk]  =rti^Tn-  y<k+1* 

*[Afc(*>]  -  ^CTk^  *1  -  ^  Et[Tk] 

*  t^TT  [kc  Y(k*  xt)  '  r  Y{k+1'  xt)] 


The  recurrence  relation: 

Y(k+l,x)  =  k  Y(k, x)  -  xk  e”X  , 

may  be  used  to  obtain  the  following  alternative  form  for  E^Ak(t)J  : 

f[Ak(t)]  =  jk(c  -  J)  Y(k,Xt) 

The  unconditional  payment  process  has  expectation  of  the  form: 
£[A(t)]  =  L  P[N=k] 


+  b  \k_1  tk  e”Xt 


k=l 


-Xt 


*(«-£)  l  •dr7Tlf0'.>-t)P[N=k]+if-^  l  iS?fyrPtN=k] 

k=l  k=l 

The  a  priori  probabilities  P[N=k]  can  be  selected  to  represent  any 
particular  kind  of  failure  model  one  might  wish  to  study.  They  should  indicate 
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what  is  perceived  to  be  the  situation  in  whatever  actual  availability  profile 

is  being  studied.  One  likely  (though  not  necessary)property  of  P[N=k] 

might  be  that  they  should  decrease  with  increasing  k  so  that  one  failure  is 

most  likely,  two  failures  less  likely,  and  so  forth.  Since  the  consideration 

is  for  a  finite  time,  the  warranty  period,  such  a  priori  considerations  seem 

reasonable.  A  specific  form  for  such  a  priori  probabilities  is: 

-k 

P[N=k]  =  (e-1)  e  Using  this  particular  example  gives: 

•  k 

I  $Mfyr  pCN=k^  = 

k=l 

so  that,  for  this  particular  example; 

•  -k 

E[A(t)]  =(c-|)(e-l)  I  Y(k,U)  +  b(l-e_1)te‘Xt  exp(U/e) 

k=l 


The  results  obtained  above  can  be  given  in  an  alternative  form  using  the 
relation: 

Y(k,x)  =  (k- 1 ) !  [l-e"x  ek4(x)] 

where  e^  k(x)  is  the  truncated  exponential  series  defined  as: 

k- 1 

,  .  V  xm 

®k-l(x)  "  L  ~  * 

m=0 


The  expected  payment  process  for  general  P[N=k]  has  the  form: 

OB 

ECA(t)]  =  (c  -£)  l  k[l-e'Xt  ek  l(Xt)]  P[N=k] 
k=  1 


+  I  ^rrrp[Nak] 


k=  1 
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E[A(t)]  =  (c  -J;  E[N]  +  e'U 


P[N=k]  . 


Section  4.  4  Limit  Results  and  Checks 

Three  kinds  of  limit  situations  may  be  considered  for  the  incentive 
warranty  models  presented  in  this  chapter.  Such  limits  provide  checks  on 
the  analysis  and  in  some  cases  give  useful  bounding  forms  for  the  results. 
They  also  allow  a  comparison  with  intuitive  impressions  about  the  model 
behavior  in  special  situations.  One  kind  of  limit  is  to  reduce  the  time 
between  the  last  two  failures  to  aero.  This  was  used  as  a  check  on  calcula¬ 
tions  in  Chapter  3  and  serves  the  same  purpose  here.  Another  type  of  limit 
is  to  remove  the  finite  time  restriction  so  that  full  expectations  can  be  used. 
Such  limit  results  give  simple  forms  that  can  be  more  easily  appreciated 
from  an  intuitive  viewpoint  than  the  complicated,  finite  time  results.  A 
third  type  of  limit  is  to  allow  all  the  failure  rates  to  become  the  same  in  the 
general  model  to  obtain  the  equal  failure  rates  case.  This  seems  to  be  a 
particularly  difficult  limit  process  and  only  partial  results  were  obtained  as 
given  in  Section  4. 4.  3. 

4.4.  1  Coalescence  of  Last  Failures 
As  an  example  of  the  coalescence  of  last  failures  consider  the  density 
function  f^(t)  which  should  become  in  the  limit  as  X^  -*•.  The  calcula¬ 

tion  follows: 
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The  last  limit  is  zero  and  the  other  limit  is  -1  so  that  result  is: 


The  general  coalescence  result  can  be  established  from  the  integral 
difference  equation: 


r  -X.  (t-w) 

£k(t)  =  '  fk-l(w)  Xk  *  dw 

0 


t 

r  i<*> 
0 


-\|^t  r  -Xj.(t-w)  1 

lim  fk(t)  =  fk  j(t)  +  lim  -e  ffc  j(0)  -j  ffc  j(w)  e  dwl  . 

V  '  Xk^"  0  J 

From  which  the  general  result  follows: 

lim  fk(t)  =  fj^  jft)  . 

V" 

The  same  argument  applies  to  the  equal  failure  rate  case  in  the  sense  that  one 
X  is  allowed  to  become  unbounded  while  all  the  others  remain  fixed  in  value. 

4.  4,  2  Unbounded  Failure  Time 

The  quantity  P[Tk  <  t]  should  of  course  approach  one  as  t  becomes 

unbounded.  It  is  the  distribution  function  F—  (t).  For  the  equal  failure  rate 

case  the  form  of  F—  shows  at  once  that  lim  F—  (t)  *  1.  This  result  is  well 
Tk  t-w  Tk 
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known  as  a  consequence  of  the  fact  that  has  a  gamma  distribution  in  the 
equal  failure  rate  case. 

The  unequal  failure  rate  expression  is: 


lim  Ft  (t)  =  \ 
t-'®  K 


k 

v  ■  i 

“W1? 


=  A, 


and  it  is  by  no  means  obvious  that  this  expression  is  equal  to  one.  Direct 
algebraic  procedures  show  that  for  k=l,2,  3  the  quantity  is  indeed  one. 
However  such  procedures  become  extremely  involved  for  larger  values  of  k. 
One  approach  is  to  show  that  the  expression  is  independent  of  the  Xj 
quantities  and  then  evaluate  it  for  a  particular  set  of  values  to  show  that  it 
is  the  constant  one.  However  it  seems  to  be  difficult  to  show  for  example 
that  the  derivative  with  respect  to  a  particular  is  zero,  from  which  the 
constant  nature  of  the  quantity  Afe  would  follow. 

Another  approach  can  be  taken  which  shows  that  A^  is  constant. 

The  quantity  A^  can  be  shown  to  be  constant  by  using  the  integral  difference 
equation  from  page  4.  8.  By  direct  verification  it  has  been  established  that 
k  -X.  t 

V* >  *  •  •  •  *i  l  a,-*.).1..  ,atis£ie‘ the  “»"atio"-  “ ,hi* form 

i_  j 

is  used  to  express  f^  j(w)  in  the  integrand  of  the  equation,  upon  integration 
f^(t)  must  result.  The  calculation  of  the  integral  follows: 


k-1 

r  *k 

c-i  L  (xf 

i=i  1 


-Xiw  -Xk(t-w) 


xi>*'#  ^k-l'tyj 


dw 


Col  ~Xr  tk-l 

V  e  Xk*  V  1 

L  <x  L  -  xt)  *  *  •  -  x^  *  e  (x1-xi)-**(xk- 
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This  must  equal  f^t)  which  can  be  expanded  as  follows: 


fk(t>  =  xr"xk 


■tj  *X1  "Xi*  ‘  “  ^Xk  “  Xi^  ^l^P 


(Xk-l  "Xk* 


Setting  the  two  expansions  for  f^(t)  equal  yields: 


(Xl'Xi)*'‘  (Xk'Xi}  ’  <X1-Xk^  (Xk-l‘Xl 


Now  consider: 


k-1 

*  V  1 

Ak-i -xr-*xk-i  l  xiix1-x.y---(xk  rx.) 


Vxi 

Ak-1  =xr-’xk-i  l  x^;-xJ"Wxk-y 


Ak-l=xl**-xk  L  \A\i-\.)-' 
i=l  1  1  1 


,  ,  V  1 _ 

-xr**xk-i  l  (Xi-x^-fXj.-Xi) 

i=l  i  1  K  1 

Using  the  result  obtained  above  from  the  integral  difference  equation  for  the 
last  summation  yields: 


Afc-l  =  xi***xk  I  Xi(Xi  -Xj)*  *  •  (Xi -XiJ  + 


Xl*  *  ’ Xk-1 


Xj(X  #  •  (Xj^X^)  (X ^Xj^)*  •  •  (Xfc.  i"X|^7  Xj^ 


Multiplication  of  numerator  and  denominator  of  the  last  term  by  Xk  as  shown 


shows  that: 
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Ak-i =  xr 


k 

r  _ i _ 

‘Xk  L 

i=  1 


=  A, 


Thus  A^  does  not  depend  on  k.  It  is  a  constant  as  required.  The  constant 
expression  can  be  evaluated  by  assigning  the  values  =  i  for  i  =  1,  •  •  * ,  k. 
The  result  is: 


' — 1  1 

Ak  =  k!  1  i(l-i)(2-iV  •  •  (k-i} 


i=l 


k!  / 


L  TITT 

i=l 


rrrmvr^ittr 


k,  v 

L-  i!  (k-i) ! 
i=l 


=  L  (-D 
i=l 


i-1  'k^ 

\  V 


IV 

Since  \  (j)(-l)1  =  °»  it  follows  that  Afc-  1=0  or  Afc=l  as  required. 


i=0 


Another  limit  consideration  for  unbounded  t  is  the  study  of  the  pay¬ 
ment  process  A(t). 

*  T  ~\  k  *T  ~! 

For  the  equal  failure  rate  case  it  was  found  that  E^T^j  =■?«•  R,  T^+j <tj 

'^T,<tJ  =  1  for  all  k.  Therefore: 


so  that  lim  E 

t^CD 


t[Tk]  *T  •inc*  lto  p| 


4.1? 


lim  E[Ak(t)  ;  =  kc  -  ^  ,  and 


lim  E[A(t)]  =  ^  lim  ELAk(t)J  P[N=k] 


=  (c  -  r)  1  k  p[N=k] 


■  (<=  -  t)  eCn]  • 

This  result  is  intuitively  satisfying,  indicating  in  the  full  expectation 
case  that  the  payment  under  the  warranty  is  expected  to  equal  the  product  of 
the  payment  per  failure  times  the  expected  number  of  failures. 

Two  special  examples  of  this  result  are  obtained  from  expressions 
derived  in  Section  4.  3.  One  form  of  the  expectation  was  found  to  be: 


e'CAW]  =  (c-£)e[N]  +  I  -  k(c-£) 

u- 1  -  J 


P[N=k] 


and  this  yields  the  above  limit  as  t-os  from  a  different  point  of  view  (i.  e.  a 
detailed  formula  for  E(t))  . 

When  P[N=k]  =  (e-1)  e  it  was  found  that: 


E[A<t)]  =  (c-£)<e-l)  2,  (^yrY(k^t)+|(e-l)t 


-Xt(l-i) 


Since  lim  Y(k,  Xt)  =  F(k)  =  (k-1)! ,  and  the  last  term  goes  to  sero  as  t 

t-*® 

becomes  unbounded,  one  obtains: 
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lim 

t—OO 


E[A(t)]  =  (c  (e-1)  L  ke 

k=l 


In  this  expression  E[N]  =  (e-1)  ^  k  e‘k  =  ss  1.  6.  It  may  be  noted  that 

k=l 

the  value  E[N]  =  1.6  for  this  special  a  priori  distribution  of  the  number  of 
failures  indicates  that  it  is  a  reasonable  one  to  assume  if  one  and  two  failures 
seem  about  equally  likely  and  more  failures  seem  unlikely. 

4,4,  3  Approaches  to  Equilisation  of  Failure  Rates 
It  would  be  interesting  to  show  that  the  limit  of  general  expressions 
for  the  unequal  failure  rates  case  tend  to  corresponding  expressions  in  the 
equal  failure  rate  case  as  the  failure  rates  approach  the  same  value.  This 
seems  to  be  a  difficult  limit  process  to  evaluate  and  only  a  little  actual  pro¬ 
gress  has  been  made  toward  its  study.  Though  several  attempts  at  taking 
the  limit  seemed  promising  none  succeeded  in  yielding  results.  The  problem 
with  the  general  analysis  is  the  complex  interaction  of  indeterminate  forms. 
For  their  study  it  is  necessary  to  deal  with  products  containing  many  factors, 
differentiation  of  which  is  extremely  involved.  On  the  other  hand  algebraic 
procedures  also  become  extremely  complex.  One  result,  forming  a  partial 
contribution  to  the  analysis,  is  the  following: 

lim  f3  (t)  =  f3(t) 

1 1*  X2*  Xj-’X 


This  result  may  be  shown  by  direct  calculations,  starting  from  the  expres¬ 
sion  for  f3(t)  in  the  form: 


*  X j  X2  Xj 


-X^  -\3t  -x2t 

e  ,  e  e 

{Xj^XjMXj-X^)  (X2  “X.  3  )(X  1^X3)  (X  j  ~X2  HX^Xj) 


Let  Xj  =  X  to  obtain: 


lim  f3(t)  =  X  X2  \3 
X 


-xt 


(X-X2)(X-X3) 


-(X  3-X)t  -(X2-X)t 

e  e 

(X2-X3)(X-X3)  "  (X-X2)(X2-X3) 


When  X2~»X  it  is  necessary  to  evaluate  the  expression: 

-(X2-X)t 
L2^X2_X3J 


A  =  lim 

X2"X 


1  e 

(X-X,)(X-X,)  '  (X-X,)(X,-X,) 


using  some  algebra  and  L’Hospital's  rule  yields: 

(l  +  (X-X3)tJ 


A  =  - 


(X-X- 


,  so  that 


r  -a3-x)t 

lim  f  <t)  =  X2X3  e‘Xt  -£ - 2 - 

X1.X2-X  L(X-X3)  (X-X3)  u 

Similar  treatment  of  this  expression  yields: 

-(X3-X)t 

lim  f-(t)  =  X3  e'Xti  lim  - 

X  i»  X2*  X3"*X,  X3~X  K~K3 


-  1 


and  finally: 


3  .It  r  - 

lim  f3{t)  =  X  e  -y  -  f,(t)  as  required. 

X1»X2,X3“*X 


One  might  hope  that  this  kind  of  procedure  could  be  extended  to  the 
general  case.  However  attempts  to  do  so  have  remained  ineffective  so  far 
in  the  present  research.  Though  interesting, such  considerations  are  not 
part  of  the  main  analysis  and  are  not  necessary  for  the  continued  develop¬ 
ment  of  the  study  of  incentive  warranty  models. 
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Complement  to  Chapter  4 

The  varification  that  f^(t)  satisfies  the  integral  difference  equation 
may  be  carried  out  as  follows:  Let 

p*  -\k(t-w) 

R  =  !  fk_j{w)  Xk  e  dw  where  it  is  assumed  that 

0 


k-1 


fk-l(w)  =  X1  ' '  *  Xk-1  L 


-X.w 

1 


i=  1 


(XA  -  Xj) 


(Xk-l  '  V 


and  X,  t  X-  for  i  =  1,  •  •  * ,  k-1. 

k  l 

Then  by  integration: 


rk- 1  -X*  t  X  t  k”  ^ 

T  e  1  V  X  _ 1 

-e  .A  rxl-xi>- ■ 


,Xk‘XiJ 


Thus  it  will  be  established  that  R  =  f^(t)  if  it  can  be  shown  that: 


k-1 


l 

i=l 


1 _  -1  _ 

(Xl‘Xi)  ’  ’  (Xk"Xi)  ^X 1  ”Xk^ ’  *  *  (Xk-l'Xk^ 


Consider  the  partial  fraction  expansion  of  the  quantity  on  the  right  to  yield 

an  algebraic  identity,  valid  for  all  values  of  Xy  (except  of  course,  that 

*  V .  for  i=l, . . . ,  k-1).  The  expansion  is 
*  1 


•  ”  -'k)-  •  •  (Xk.1-Xk)  (XL-Xk)  (xk-l'xk) 


•rr  arator  of  both  sides  must  be  equal  so  that: 

^  :  '*2*Xk) '  ’  ’  (Xk-l"Xk>  +  B2(XrXk) '  *  '  (Xk-l”Xk* 

*  ...  *  Bk_1(X1-Xk)  •  ••  (Xk.2-Xk) 
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Let  to  obtain: 


-1 


1  (X2 ”X  j )(X  j “X  j.)  •  •  •  (Xj^_  ^ — X 


Let  X^  =  X2  to  obtain: 


B,  = 


-1 


2  '  (X1-X2)(X3-X2)***  ^k.i‘^2)  * 
And  so  forth,  finally  letting  X^  =  X^  2  to  obtain 

B, 


-1 


k_1  <Xl"Xk-l)***  (Xk-2"Xk-l) 


Putting  in  these  values  and  interchanging  the  last  factor  in  each  denominator, 
thereby  changing  the  sign,  yields  the  desired  summation.  This  establishes 
f^(t)  as  the  solution  form  for  the  integral  difference  equation  (page  4.  8). 
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Chapter  5  Extensions 

The  study  of  incentivised  availability  warranties  carried  out  as 
described  in  this  report  suggest  a  number  of  additional  investigations.  These 
range  from  technical  and  theoretical  questions  to  implementation  of  the  meth¬ 
odology  for  structuring  warranty  contracts.  Various  extensions  of  the  method- 
logical  framework  may  also  be  identified.  A  number  of  these  considerations 
are  set  forth  in  this  chapter  as  a  matter  of  record  and  to  serve  as  possible 
guidelines  for  subsequent  research. 

Applicability 

The  very  concept  of  incentive  availability  warranties  implies  application 
to  actual  contract  structure  either  directly  or  in  some  modified  form.  Its 
analysis  has  no  meaning  except  as  related,  in  some  context,  to  application. 

The  present  theoretical  introduction  of  a  methodology  for  study  of  these  con¬ 
cepts  may  be  carried  toward  application  in  several  ways.  The  rather  complex 
forms  obtained  for  variance  and  expectation  may  be  investigated  in  numerical 
(expressed  as  graphical)  terms.  The  variance  is  a  measure  of  risk  assumption 
and  as  such  plays  a  central  role  in  sophisticated  contract  structuring.  However 
the  formulas  for  variance  in  the  important  one  and  two  failure  situations  are 
complicated  functions  of  the  contract  parameters,  difficult  to  appreciate  in 
terms  of  analytic  expression.  Numerical  investigation  of  these  formulas  and 
of  the  expected  values  for  both  special  and  general  situation  formulas  will 
carry  the  methodology  further  toward  actual  representational  analysis  of  incen¬ 
tive  warranty  structures.  A  major  goal  of  contract  structure  analysis  is  to 
study  the  interplay  between  motivation  of  contractor  toward  customer  goals  and 
the  division  of  risk  assumption  in  a  potentially  profitable  environment.  Further 
detailed  study  of  the  methodological  formulas  will  contribute  to  this  basic  con- 
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side  ration. 


Another  feature  related  to  application  of  the  theory  to  actual  warranty 
development  is  the  interplay  between  the  kinds  of  risks  a  contractor  must 
asstime.  In  a  design  decision  context  the  contractor  may  expend  money  in  the 
manufacturing  phase  to  reduce  the  risk  of  spending  money  in  the  warranty 
period.  Incentivisation  of  warranties  would  reduce  the  motivation  for  a  con¬ 
tractor  to  so  reduce  warranty  cost  and  hence  work  counter  to  high  availability. 
On  the  other  hand  the  incentivisation  may  motivate  the  contractor  to  agree  to 
a  warranty  contract.  Thus  the  situation  is  complex  for  both  parties  and  would 
doubtless  benefit  from  further  investigation. 

Applicability  studies  are  therefore  of  two  distinct  types:  investigation 
of  how  desirable  the  incentivising  concept  actually  is  in  various  contexts,  and 
detailed,  numerical  study  to  provide  contextual  appreciation  of  the  methodo¬ 
logical  results. 

Extensions  of  Methodology 

There  are  a  number  of  considerations  that  would  extend  the  potential 
application  of  the  methodology  by  making  it  representational  of  a  broader  class 
of  situations.  Some  of  these  deal  with  the  incentive  pool  structure,  some  with 
the  underlying  stochastic  environment,  and  some  with  the  model  structure  itself. 
A  number  of  these  extensions  are  briefly  identified  below: 

•  Change  in  the  form  of  the  incentive  fee  pool  may  allow  alternative 
interpretations  and  anlayses.  As  a  simple  illustration  an  exponential  type  fee 
pool  may  be  used.  In  the  single  failure  (reliability)  case  the  payment  process 
has  the  form  A^(t)  =  c  -  aQ  e^*.  With  exponential  failure  this  model  gives: 

dividing  nicely  into  the  payment  without  incentive  pool  and  the  contribution  made 
by  the  incentive  fee. 


r 


•  The  fee  pool  form  can  change  after  a  failure  rather  than  always 
return  to  the  same  rate  of  increase.  The  change  can  a  fixed  value,  depending 
on  the  failure  time,  or  be  a  random  variable  the  value  of  which  is  specified 
by  a  probability  distribution  assumed  as  part  of  the  model. 

•  Alternative  contract  forms  can  be  used  to  produce  motivation 
similar  to  the  incentive  contract  approach.  For  example  the  award  fee  con¬ 
cept  can  be  applied  to  warranty  contracting.  Awards  can  be  designed  to  off¬ 
set  contractor  payments  under  the  warranty  (as  the  incentive  pool  does). 

They  can  be  based  on  the  occurrence  of  failure  events,  the  nature  of  the 
failure,  or  the  time  between  observed  failures. 

•  In  the  exponential  failure  model  some  of  the  failure  rates  can  be  the 
same  and  others  different.  This  is  a  potentially  useful  model  falling  between 
the  two  extremes  considered  in  Chapter  4.  In  particular  it  would  be  interesting 
to  study  one  case,  e.g.  for  k=3,  in  detail.  Such  a  study  may  be  expected  to 
shed  insight  on  the  general  case. 

•  The  failure  rates  can  be  considered  as  variable,  leading  to  more 
general  failure  distributions  such  as  the  Weibull. 

•  In  the  present  study  uniform  failure  time  distributions  were  used  only 
for  the  one  failure  case.  However  these  distributions  provide  detailed  instru¬ 
ments  for  some  investigations  of  the  payment  process.  Their  use  can  be 
extended  to  more  than  one  failure  situations. 

•  The  cost  to  repair /replace  an  item  need  not  be  constant.  Variable 
service  cost  can  be  of  different  kinds,  either  being  selected  from  a  (small) 
number  of  discrete  values  or  as  a  value  of  a  continuous  random  variable.  An 
interesting  situation  is  when  cost  is  represented  by  a  mixed  random  variable 

in  which  a  fixed  replacement  cost  has  non  sero  probability  and  a  range  of  repair 
costs  are  continuously  distributed.  Consideration  of  replacement/ repair  time 


and  cost  will  increase  the  potential  applicability  of  the  methodology  and 
indicate  the  effect  of  such  factors  on  the  concepts  of  incentivised  warranties. 

•  The  time  to  repair /replace  can  be  explicitely  considered  in  the 
methodology.  The  service  time  can  be  thought  of  as  a  single  random  variable 
(deterministic  service  time  is  a  special  case)  or  for  more  detailed  study  it 
can  be  represented  as  a  sum  of  waiting  time,  diagnosis  time,  and  actual  repair 
(replacement)  time.  Its  consideration  alone  can  represent  an  extensive  area 
for  research  and  it  is  necessary  to  investigate  only  the  general  effect  of 
service  time  on  the  availability  profile  as  it  relates  to  the  payment  process. 

•  The  concepts  of  availability  warranties  can  be  extended  to  logistic 
warranties  in  which  additional  features  of  the  logistic  process  are  brought 
into  the  methodological  considerations.  An  instance  of  such  additions  is  the 
multiple  equipment  warranty  in  which  multiple  failures  can  occur  and  require 
service.  The  service  queue  enters  into  the  considerations  as  does  the  full 
logistic  pipeline.  Analysis  involving  such  ideas  should  be  limited  to  con¬ 
siderations  having  the  strongest  impact  on  the  payment  process.  It  seems 
likely  that  features  of  the  service  queue  (such  as  the  number  of  service 
channels  and  the  queue  discipline)  may  be  important  in  this  context. 

•  Considerations  related  to  the  value  of  money  are  important  in 
applying  the  incentive  warranty  methodology  to  actual  contracts.  These  con¬ 
siderations  are  distinct  from  the  availability  theory  approach  to  the  model 
methodology  but  are  associated  with  it  in  the  formulation  of  realistic  contracts. 
There  are  many  distinct  possibilities  for  dealing  with  incentive  fee  pool  dollars 
and  contractor  payments  from  the  discounting  point  of  view. 

Discounting  can  simply  be  incorporated  into  the  model  as  a  correction 
calculation  on  the  value  of  money. 

A  more  interesting  situation  is  when  the  contractor  can  select  less  than 
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the  available  fee  pool  money  and  retain  some  toward  the  future  failure  expense. 
When  the  value  of  money  is  decreasing  this  seems  like  a  bad  idea  since  the 
full  amount  available  will  never  be  worth  more  to  the  contractor  than  when  it 
becomes  available.  On  the  other  hand  the  contractor  may  have  an  alternate 
use  for  the  money  that  may  increase  its  value  in  the  future.  It  can  be  seen  that 
there  are  many  possibilities.  To  be  used,  they  must  be  allowed  under  contract 
specification  and  their  motivational  effect  on  the  contractor  must  be  carefully 
studied.  If  conducted  in  detail  such  considerations  become  an  additional  dimen¬ 
sion  for  the  incentivised  warranty  methodology. 

Technical  Considerations 

The  result  of  allowing  different  failure  rates  to  become  equal  is  a 
particularly  difficult  limit  evaluation  that  would  contribute  to  the  internal  con¬ 
sistence  of  the  research. 

These  things  are  not  required  for  the  main  line  of  methodological 
development  but  do  present  some  technical  features  of  the  mathematical  analy¬ 
sis  that  could  benefit  from  further  study. 

Theoretical  Questions 

The  use  of  what  has  been  called  partial  expectations  in  this  report 
introduces  questions  about  the  role  of  this  concept  within  the  general  frame¬ 
work  of  applied  probability.  It  is  not  a  conditional  expectation  though  it  is 
used  to  compute  conditional  expectations  in  the  usual  sense,  e.g.  EtAk<*»] 
a  conditional  expectation  whereas  Ej  T^j  is  not.  Some  study  of  this  concept 
has  been  given  in  Section  4.  1.2. 

The  other  major  theoretical  consideration  is  the  use  of  a  priori 
probabilities  for  P[N=k]  and  the  relation  between  such  probabilities  and 
those  generated  by  the  failure  time  distributions.  Because  of  the  time  bound¬ 
ary  imposed  by  the  warranty  period  the  failure  process  does  not  develop  under 
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the  operation  of  the  underlying  stochastic  process.  The  point  of  view  is  rather 
one  in  which  the  number  of  failures  is  specified  and  act  to  condition  the  process 
leading  to  the  use  of  partial  expectations  of  failure  time  and  conditional  expec¬ 
tations  of  payment  up  to  a  time  t  less  than  the  warranty  period.  However  it 
may  be  that  the  usual  process  for  generating  failures  should  play  a  role  in  the 
incentive  availability  warranty  methodology. 

In  the  general  availability  model  the  number  of  failures  generated  by 
the  stochastic  failure  process,  Nfc,  is  a  random  variable  of  interest  in  addition 
to  the  failure  time  distribution.  This  random  variable,  generated  by  the 
process,  is  distinct  from  N  the  number  of  failures  postulated  for  the  equipment 
as  an  a  priori  feature  of  the  model.  The  two  variables  Nfc  and  T,  one  discrete 
and  one  continuous,  complement  each  other  as  representational  factors  of 
the  random  process.  If  T  and  Nt  are  the  failure  time  and  number  of  failures 
to  time  t  respectively  the  relation  between  them  is  expressed  by  the  equation 
FT(t)  =  1  -  p("Nt  =  Oj.  By  differential  difference  equations  the  probability  dis¬ 
tribution  for  Nt  can  be  determined  from  F^,(t),the  distribution  function  for  T  . 
In  such  relations  some  assumptions  about  the  process  must  also  be  made  (such 
as  independence,  or  non-cascade  assumptions).  The  study  of  the  role  of 
contributes  to  the  understanding  of  any  stochastic  process.  This  is  true  of  the 
payment  process  model  of  incentive  availability  warranties,  where  there  is  the 
added  feature  of  possible  application  of  the  model  to  actual  contract  structuring. 
Potential  application  is  enhanced  by  being  able  to  consider  the  model  from 
different  viewpoints. 
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